A0-A6r?  5t»  ENVIRONMENTAL  RESEARCH  INST  OF  MICHIGAN  ANN  ARBOR  IN— ETC  F/G  17/5 
STATISTICAL  ANALYSIS  OF  TERRAIN  BACKGROUND  MEASUREMENTS  DATA*  <U> 

MAR  77  R  SPELLICY  ,  J  BEARD  »  J  R  MAXWELL  N00123-76-C-0708 

UNCLASSIFIED  ERIM-120500-12-F  NL 


;  FILEJQEL 


I  Final  Report 

STATISTICAL  ANALYSIS  OF  TERRAIN 
BACKGROUND  MEASUREMENTS  DATA 


infrared  and  Optics  Division 


MARCH  1977 


Approved  (or  Public  Release : 
Distribution  Unlimited 


Optical  Signatures  Program 
Naval  Weapons  Center 
China  Lake,  California 


Notices 


Sponsorship  The  work  reported  herein  was  conducted  by  the  Environmental 
Research  Institute  of  Michigan  (formerly  the  Willow  Run  Laboratory  of 
The  University  of  Michigan)  for  the  Naval  Weapons  Center,  China  Lake, 
California  under  Contract  Number  N00123-76-C-0708.  The  Project  Managers 
were  Dr.  Lowell  Wilkins  and  Mr.  Ted  Smith. 

Disclaimer  The  views  and  conclusions  contained  in  this  document  are 
those  of  the  authors  and  should  not  be  interpreted  as  necessarily 
representing  the  official  policies,  either  expressed  or  implied,  of 
the  Department  of  Defense  or  of  the  U.  S.  Government. 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


UNCLASSIFIED 


PERFORMING  ORGANIZATION  NAME  »NO  ADDRESS 

Infrared  and  Optics  Division  v-^ 

Environmental  Research  Institute  of  Michigan 
P  0  Box  618,  Ann  Arbor,  MI  48107  _ _ 

CONTROLLING  OFFICE  NAME  ANO  ADORES*  /\  / 

Dr.  Lowell  Wilkins,  Code  51403  (713  939-3814) ' Vi- 
Naval  Weapons  Center  _ -- 


IO  PROGRAM  ELEMENT  PROJECT  TASK 
AREA  4.WSIW  UNIT  NUMBERS 


14.  MONITORING  AGENCY  NAME  ANO  AOORESf 
(if  different  from  Controlling  Off  tit ) 

Receiving  Officer 
Naval  Weapons  Center 
China  Lake,  CA  93555 


IB  DISTRIBUTION  STATEMENT  (of  Ikh  R  tfnrt) 


13  NUMBER  OF  PAGES 

307 


IS  SECURITY  CLASS  (of  tmn  npn) 

Unclassified 

13a  DECLASSIFICATION/ OOWNGRAOl 
SCHEOULE 


Approved  for  Public  Release:  Distribution  Unlimited 


17  DISTRIBUTION  ST  ATE  ME  NT  (of  the  aburoct  mitred  in  Block  20.  if  different  from  Report) 


,  on/mue  on  rerene  ode  if  net  earn  ry  and  identify  by  block  number) 

Terrain  backgrounds  data,  statistical  analysis,  multispectral  data 


A8S7 XjACT  (Continue  on  rerene  ude  if  *  ternary  and  idjwifj  by  block  number) 

Several  terrain  backgrounds  scenes  were  selected  and  statistics  derived 
from  calibrated  multispectral  airborne  scanner  data.  Conventional  statis¬ 
tical  parameters  -  means,  standard  deviations,  histograms,  spectral  cor¬ 
relations,  and  power  spectra  -  are  reported  for  the  1 . 0-1 . 4  jjif,  1.5-1. 8  (aaP;~ 
2.0  -  2.6 4. 5-5. 5  yrfj  and  9.3-11. 7  spectral  bands.  In  addition, 
a  program  has  been  Initiated  to  develop  equivalent  elliptical  area/intensity 
statistics  for  these  background  scenes  which  more  adequately  describe  — »>  ^ 


EDITION  OF  I  NOV  *3  IB  OBSOLETE 


UNCLASSIFIED 

SECURITY  classification  OF  THIS  PAGE 


/Kic  (  £>  vx* 


v  ^  i 


UNCLASSIFIED 


EC^^T' 


security  CLASSIFICATION  OF  THIS  PAGE  (U'kn  D*i*C*iirtJ) 


the  infrequent  occurrence  of  small  areas  of  high  intensity 
seekers  and  trackers  false  alarm  problems.  Area/ intensity 
are  reported  for  two  of  the  scenes. 


\ 


which  cause 
statistics 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  FAGE  fVkm  0*1.1  tlltmi) 


M>C  TAB 
Ifcaomom 
Just  if i  . 


TABLE  OF  CONTENTS 

Introduction  and  Sunmary  .  . . 7 

The  Multispectral  Scanner  ...  .  H 

Procedures  and  Selected  Examples . 23 

3.1  Preprocessing  of  Scanner  Data . 23 

3.2  Statistics  Definitions . 27 

3.2.1  Point  Statistics:  Means  and  Standard  Deviations  .  .27 


3.2.2  Correlative  and  Area  Statistics 


3.3  Modeling  the  Background  Scene . 30 

4.  Discussion  of  Results . 57 

5.  Summary  and  Recommendations  . 61 

Appendix  I:  Details  of  Calibration  and  Field-of-View 

Averaging  Procedures  .  1-1 

Appendix  II:  Complete  Compilation  of  Statistical  Data  . II-2 


LIST  OF  FIGURES 


2.1  Airborne  Multispectral  Scanner  Operation  .  12 

2.2  Optical  Schematic  of  ERIM  Experimental  Multispectral 

Scanner,  M-7 .  14 

2.3  ERIM  Experimental  Multispectral  Scanner  System  .  16 

2.4  Scanner  Voltage  Output  Versus  Time  .  18 

2.5  ERIM  M-5  Multispectral  Scanner .  21 

3.1  Schematic  Representation  of  the  Procedures  used  for 

Field-of-View  Averaging  .  26 

3.2  Subarea  Histograms  for  the  9.3  -  11.7  pm  Channel  of  Flint-1  .  30 

3.3  Total  Area  histogram  for  the  9.3  -  11.7  pm  Channel  of  Flint-1  31 

3.4  Equivalent  Elliptical  Areas  for  the  9.3  -  11.7  pm  Channel  of 

Flint-1  with  a  1.5a  Threshold .  36 

3.5  Equivalent  Elliptical  Areas  for  the  9.3  -  11.7  pm  Channel  of 

Flint-1  with  a  2.0a  Threshold . 37 

3.6  Equivalent  Elliptical  Areas  for  the  9.3  -  11.7  pm  Channel  of 

Flint-1  with  a  2.5a  Threshold .  38 

3.7  Equivalent  Elliptical  Areas  for  the  9.3  -  11.7  pm  Channel  of 

Flint-1  with  a  3.0a  Threshold .  39 

3.8  In-Track  One  Dimensional  Averaged  Wiener  Spectrum  for  the 

9.3  -  11.7  pm  Channel  of  Flint-1 .  47 

3.9  Cross-Track  One-Dimensional  Averaged  Wiener  Spectrum  for 

the  9.3  -  11.7  pm  Channel  of  Flint-1 .  48 

3.10  Two-Dimensional  Wiener  Spectrum  for  the  9.3  -  11.7  pm  Channel 

of  Flint-1 .  49 

3.11  Artificial  Color  Image  of  the  9.3  -  11.7  pm  Channel  of  Flint-1  51 

3.12  Artificial  Color  Representation  of  the  Pseudo-Image  Generated 

From  the  9.3  -  11.7  pm  Channel  of  Flint-1 .  51 


4 


LIST  OF  TABLES 


2.1  M7  Scanner  Performance  Characteristics  .  ....  19 

2.2  Detector  Configurations  for  ERIM  M7  Scanner .  20 

3.1  Sub-Area  and  Total  Area  Means  and  Standards  Deviations 

for  Flint-1 .  29 

3.2  Correlation  Coefficients  for  the  Reflective  IR  Channels 

in  Sub-Area  2  of  Flint-1 .  34 

3.3  Area  Distribution  for  Recognized  Regions  in  the  9.3  -  11.7  pm 

Channel  of  Flint-1  with  a  lo  Threshold .  41 

3.4  Perimeter  Distributions  for  Recognized  Regions  in  the 

9.3  -  11.7  pm  Channel  of  Flint-1  with  la  Threshold  .  42 

3.5  Shape  Factor  Distributions  for  Recognized  Regions  in 

the  9.3  -  11.7  pm  Channel  of  Flint-1  with  a  lo  Threshold  .  .  43 

3.6  Area  Distributions  for  Recognized  Regions  in  the  Real 

and  Pseudo-Image  of  the  9.3  -  11.7  pm  Channel  of  Flint-1  .  .  53 

3.7  Area  Distributions  for  Recognized  Regions  in  the  Real 

and  Pseudo-Images  of  the  9.3  -  11.7  pm  Channels  of  Flint-1  .  54 

: 

3.8  Area  Distributions  for  Recognized  Regions  in  the  Real 

and  Pseudo-Images  of  the  9.3  -  11.7  pm  Channels  of  Flint-1  .  55 

3.9  Area  Distributions  for  Recognized  Regions  in  the  Real 

and  Pseudo-Images  of  the  9.3  -  11.7  pm  Channels  of  Flint-1  .  56 

4.1  Total  Area  Data  Summary . . .  58-59 


5 


PRECEDING  PAG*  BLAJHUNOT  FILMED 


l 

1 

INTRODUCTION  AND  SUMMARY 

There  is  a  need  for  both  target  and  background  data  for  the  design 
of  sensors  to  detect  and  track  targets  against  sky,  cloud,  and  terrain 
backgrounds.  The  objective  of  this  backgrounds^analysis  program  is  to 
develop  backgrounds  data  useful  to  systems  designers.  The  first  phase 
of  the  program,  which  has  been  completed  [Reference  1],  was  to  survey 
the  available  backgrounds  data  and  compile  an  index  -  bibliography  of 
pertinent  information.  In  general  there  is  considerable  data  on  ter¬ 
rain  backgrounds  but  much  less  on  sky  and  cloud  backgrounds.  Data  is 
especially  lacking  for  two-dimensional  high  spatial  resolution  data 
on  clouds  in  the  spectral  bands  of  interest.  In  this  second  phase  of 
the  program  we  have  selected  a  variety  of  terrain  backgrounds  for  which 
high  spatial  resolution  multispectral  data  are  available.  A  variety  of 
statistical  measures  have  been  derived  from  these  data  and  are  presented 
in  this  technical  report.  These  include  the  conventional  statistical 
parameters  of  means,  standard  deviations,  histograms,  Wiener  (power) 
spectra,  and  spectral  correlations  as  well  as  new  area/ intensity  statis¬ 
tics  which  are  particularly  relevant  in  view  of  recent  advances  in  sen¬ 
sor  and  processor  technology.  Additional  efforts  to  assess  the  utility  of 
the  various  background  statistical  measures  will  be  reported  in  future 
technical  reports. 

Conventional  background  statistics  have  in  the  past  been  adequate 
since  in  many  instances  there  is  a  high  contrast  between  the  target  and 
background.  In  such  instances,  the  highest  intensity  background  points 
determine  the  highest  threshold  setting  necessary  to  eliminate  all  false 
alarms  and  the  histograms  provide  estimates  of  how  the  detection  proba¬ 
bility  and  false  alarm  rate  vary  with  threshold  setting.  However,  today 
there  is  a  need  for  higher  order  background  statistics  because  of  the 
increased  sophistication  of  background  rejection  techniques  employed 

[1]  J.  Beard,  J.  Braithwaite,  R.  Turner,  Infrared  Background  Survey  and 
Analysis,  Environmental  Research  Institute  of  Michigan,  Ann  Arbor, 

June  1976. 
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with  large  detector  arrays  and  imaging  or  scanning  sensors.  Although 
the  power  spectrum  is  a  background  statistic  which  does  vary  with  the 
spatial  distribution  of  radiances  in  the  scene,  it  is  clearly  an 
inadequate  background  descriptor  for  most  of  today's  problems.  Small 
areas  of  high  radiance  produce  most  of  the  false  alarms  in  today's 
sensors,  and  the  power  spectrum  does  not  distinguish  between  scenes 
with  many  low  intensity  areas  and  those  with  a  few  high  intensity  areas. 
This  is  a  well  known  fact  and  a  result  which  is  evident  from  the  back¬ 
ground  modeling  work  of  R.  Clark  Jones  and  the  early  working  group  on 
infrared  backgrounds,  WGIRB  [Reference  2]. 

Hence,  in  addition  to  the  conventional  statistics  on  terrain 
backgrounds,  we  have  also  developed  area/intensity  statistics  as  a 
statistical  measure  that  is  more  directly  useful  to  the  sensor  designer 
in  estimating  detection  probabilities  and  false  alarm  rates  with 
today's  sensor  and  processor  technology.  The  statistics  developed  are 
the  probabilities  that  regions  (of  various  sizes,  shapes,  and  orienta¬ 
tions)  will  occur  in  the  scene  above  a  specif ied radiance  threshold. 

The  region  descriptors  are  area,  major  and  minor  axes,  and  the  angular 
orientation  for  an  elliptical  area  that  is  equivalent,  in  geometric 
area  and  ratio  of  second  moments,  to  each  contiguous  region  above  the 
radiance  threshold.  For  most  scenes  the  number  of  regions  and  their 
areas  decrease  as  the  threshold  is  raised,  while  the  number  of  small 
regions  above  any  preset  threshold  varies  from  one  scene  to  the  next. 
These  area/intensity  statistics  are  analogous  to  the  more  familiar 
pulse  length  statistics  for  one-dimensional  records. 

The  statistical  parameters  for  the  occurrence  of  areas  and  inten¬ 
sities  in  the  various  background  scenes  are  not  only  directly  useful 
for  estimating  sensor  performance  but  may  also  be  useful  in  simulating 
whole  classes  of  backgrounds.  The  equivalent  ellipses  at  each  threshold 
can  be  positioned  to  simulate  the  actual  scene  from  which  they  were 

[2]  R.  Clark  Jones,  Dr.,  The  Study  of  IR  Backgrounds  by  the  Wiener 
Spectrum  Method,  Polaroid  Corporation,  Cambridge,  Massachusetts, 
December  1959. 
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derived,  or  repositioned  at  random  to  simulate  many  scenes  having  the 
same  area/intensity  statistics.  Such  simulations  of  backgrounds  do 
reproduce  many  of  the  spatial  characteristics  of  the  original  scene 
as  is  shown  by  example  in  this  report. 

Seven  scenes  were  selected  for  analysis  in  this  program  repre¬ 
senting  a  wide  range  of  backgrounds;  two  residential  scenes  (Flint-1 
and  Baltimore),  an  industrial  scene  (Flint-2),  a  mountainous  structured 
terrain  (Mill  Creek,  OK),  a  natural  busy  terrain  of  trees  and  hills 
with  numerous  shadows  (Black  Hills,  S.  D.)  and  two  mountainous  unstruc¬ 
tured  terrains  (Pisgah  Crater,  CA  and  Mono  Lake,  CA) .  These  were 
collected  with  the  ERIM  airborne  multispectral  scanner  system  with 
a  nominal  5  to  10  foot  spatial  resolution  in  spectral  bands  including 
1.0-1. 4  pm,  1.5-1. 8  pm,  2. 0-2. 6  pm,  4. 5-5. 5  pm,  and  9.3-11.7  pm. 
Conventional  statistics  including  means,  standard  deviations,  histograms, 
and  spectral  correlations  are  reported  here  and  calibrated  data  tapes 
are  available  to  qualified  users  for  additional  statistical  analysis 
and  simulation.  In  addition  to  these  statistics,  Wiener  (power)  spectra 
and  area/intensity  statistics  are  reported  for  two  spectral  bands  of 
Flint-1  and  Mill  Creek.  An  actual  simulation  of  the  Flint-1  scene  in 
the  9.3-11.7  pm  spectral  band  has  also  been  created  with  equivalent 
ellipses  and  is  compared  with  the  true  background  scene. 

The  remainder  of  this  report  is  broken  down  into  four  sections: 
Section  2  presents  the  ERIM  calibrated  multispectral  airborne  scanner 
system;  Section  3,  the  procedures  used  for  calibration  and  statistical 
analysis  with  selected  examples;  Section  4,  a  brief  discussion  of  the 
results  of  the  analysis;  and  Section  5,  a  summary  with  recommendations. 
The  complete  data  and  statistics  package  is  contained  in  Appendix  II. 
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THE  MULTISPECTRAL  SCANNER* 

Two  multispectral  scanner  systems  have  been  in  use  at  ERIM  since 
1968.  The  newer  M-7  scanner  was  used  at  the  time  the  Flint,  Baltimore, 
and  Mill  Creek  data  was  generated  while  its  predecessor  the  M-5  scanner 
was  used  in  gathering  the  Pisgah  Crater,  Black  Hills,  and  Mono  Lake 
data.  These  two  scanners  are  similar  so  that  only  the  M-7  will  be 
discussed  in  detail  with  the  difference  between  the  two  systems 
elaborated  on  at  the  end  of  this  section. 

The  M-7  scanner,  covering  a  wavelength  range  from  0.33  to  14.0 
micrometers,  can  operate  in  up  to  19  different  bands  of  the  ultraviolet, 
visible,  and  infrared  regions.  Of  these  bands,  12  can  be  selected  for 
tape  recording  at  any  one  time  on  a  14-track  analog  tape  machine.  More 
recently  a  digital  recording  svstem  has  been  added.  As  many  as  five 
separate  radiation  reference  sources  may  be  recorded  sequentially  along 
with  the  ground  video  once  each  scan  line.  The  total  system,  including 
boresight  cameras,  is  usually  operated  in  a  Douglas  C-47  aircraft. 

The  simplified  diagrams  of  Figure  2.1  illustrate  a  typical  line 
scanner  and  its  method  of  airborne  use.  As  shown  in  the  optical 
schematic  at  the  top  of  the  figure,  the  scanner  basically  consists  of 
an  optical  telescope  with  its  narrow  field  of  view  redirected  by  a 
rotating  flat  mirror.  This  mirror  causes  the  system  to  scan  in  a  plane 
perpendicular  to  the  longitudinal  axis  of  the  aircraft.  A  radiation 
detector  in  the  focal  plane  of  the  telescope  converts  the  focused  beam 
of  radiation  to  an  electrical  signal.  The  optical  system's  field-of- 
view  (ground  resolution  element)  first  scans  laterally  across  the  air¬ 
craft  ground  track  through  an  opening  in  the  bottom  of  the  aircraft. 

Then  before  making  the  next  ground  scan,  it  scans  radiation  references 
(not  shown)  which  are  internal  to  the  scanner.  By  the  time  the  next 
scan  begins,  the  aircraft  has  moved  forward,  thus  subsequent  line  scans 
build  upon  one  another  to  produce  a  continuous  strip  image  of  the  terrain 
beneath  the  aircraft. 

*  The  ERIM  Airborne  Multispectral  Data  Collection  is  described  in  the 
ERIM  190901-1-F  report  prepared  by  Philip  G.  Hasell,  Jr.  under  contract 
NAS  9-9304.  The  multispectral  scanner  as  described  in  that  report,  is 
presented  here  in  order  to  familiarize  the  reader  with  the  multispectral 
scanner  system  used  to  collect  the  data  which  is  being  processed  on  this 
program. 


Training  Set— 


Scan  Line 


The  multispectral  scanner  evolved  from  this  single-channel  scanner 
concept.  This  evolution  required  replacement  of  the  single  detector 
element  with  a  system  of  beamsplitters,  dispersing  optics,  and  spectral 
filters.  Figure  2.2  shows  the  optical  configuration  of  the  current 
M-7  multispectral  scanner.  A  key  feature  in  this  design  is  its  flexi¬ 
bility  for  accepting  different  radiation  reference  sources  and  new 
detector  assemblies.  Weight  and  space  savings  were  sacrificed  to 
provide  this  flexibility,  which  allows  increased  opportunities  for 
adaptation  to  a  diverse  number  of  data  gathering  modes.  Such  flexi¬ 
bility  is  an  important  attribute  for  a  general-purpose  experimental 
system. 

The  radiation  intercepted  by  the  5-inch-diameter  collecting 
aperture  is  directed  into  the  Dall-Kirkham  telescope,  which  has  a 
3-inch-diameter  secondary  mirror.  The  incoming  radiation  prevented  from 
entering  the  telescope  by  this  secondary  mirror  is  directed  upward  by 
a  folding  mirror  to  Detector  Position  1.  This  3-inch-diameter  col¬ 
lecting  aperture  operates  over  the  broad  band  of  0.3  to  14.0  pm.  To 
provide  thermal  data  at  this  position,  a  focusing  lens  designed  for  the 
8.0-14.0  pm  band  is  used  in  combination  with  a  cooled  HgCdTe  detector. 

A  dichroic  mirror  mounted  ahead  of  this  lens  diverts  ultraviolet  and  vis¬ 
ible  radiation  onto  a  photomultiplier  detector  which  is  filtered  so 
the  energy  it  receives  for  recording  is  restricted  to  a  narrow  pre¬ 
selected  band. 

The  radiation  collected  by  the  effective  4-inch  aperture  of  the 
Dall-Kirkham  telescope  is  folded  into  a  dichroic  mirror  which  reflects 
radiation  below  1.0  pm  but  transmits  that  of  longer  wavelengths.  The 
radiation  thus  transmitted  is  focused  onto  three  separately  filtered 
indium  arsenide  detector  elements  in  Position  2  by  lens  achromatized 
for  the  1.0  to  2.6  pm  region.  This  dichroic  and  lens  can  be  readily 
changed  for  different  detector  configurations. 

Radiation  at  wavelengths  shorter  than  1.0  pm  is  focused  onto  the 
entrance  slit  of  a  prism  spectrometer  at  Detector  Position  3.  The 
spectrometer  divides  and  directs  visible  and  near-infrared  radiation 
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FIGURE  2.2  OPTICAL  SCHEMATIC  OF  ERIM  EXPERIMENTAL  MULTISPECTRAL  SCANNER,  M7 


through  a  fiber-optic  image  slicer  to  as  many  as  twelve  photomultiplier 
tubes.  (In  the  current  configuration  the  radiation  goes  to  nine  separ¬ 
ate  photomultipliers.) 

The  radiation  reference  sources  are  positioned  in  line  with  the 
scan  mirror,  so  that  each  source  is  "seen"  and  registered  sequentially 
once  each  scan  line.  Currently,  five  reference  sources  are  being  used: 
an  NBS  lamp  packaged  to  simulate  a  point  source;  one  ambient  and  two 
temperature-controlled  graybody  thermal  references  that  fill  the 
collecting  aperture;  and  a  sky  illumination  reference  consisting  of  an 
opal  glass  diffusing  plate  mounted  in  the  top  of  the  aircraft.  Through 
electronic  control  of  the  lamp  and  graybodies  and  by  means  of  attenu¬ 
ating  optical  filters  for  the  sky  illumination,  the  radiation  from  all 
but  the  ambient  temperature  reference  sources  is  under  operator  control. 
During  data  collection,  all  internal  sources  are  monitored  and  recorded 
manually,  by  the  operator.  To  assure  their  validity  as  references, 
these  sources  are  calibrated  periodically  against  external  standards 
in  the  laboratory. 

The  complete  airborne  scanner  system  is  diagrammed  in  Figure  2.3. 
Terrain  radiation  enters  the  scanner  at  the  bottom  left;  radiation 
detectors  in  the  scanner  assembly  register  this  input  along  with  that 
of  the  reference  sources.  The  electrical  signals  comprising  detector 
video  outputs  are  amplified  in  preamplifiers  before  being  transmitted 
to  the  operator  console  where  the  operator  monitors  them  and  adjusts 
amplifier  gain  to  the  proper  level  for  tape  recording.  To  confirm 
satisfactory  recording,  he  is  also  able  to  monitor  signals  reproduced 
from  the  tape  record.  The  system  linearly  transforms  input  radiation 
to  voltage  analogs  which  are  recorded  on  the  magnetic  tape.  The  scanner 
system  can  generate  video  signals  in  up  to  nineteen  different  spectral 
bands  over  a  wavelength  range  extending  from  0.33  to  14.0  pm.  Any 
twelve  of  these  bands  may  be  tape  recorded  at  any  one  time  on  a  14-track 
analog  tape  machine;  the  other  two  tape  recorder  tracks  are  used  for 
housekeeping  purposes. 


The  airborne  system  (Figure  2.3)  also  includes  an  array  of  bore- 
sight  cameras  utilizing  various  film-filter  combinations.  These  aerial 
cameras  produce  film  records  often  useful  in  the  subsequent  analysis  of 
the  scanner  data. 

Electrical  voltage  representations  of  single  line  scans  for  the 
thermal  and  non-thermal  wavelength  bands  are  shown  in  Figure  2.4.  Note 
that  although  the  detectors  in  all  positions  view,  in  sequence,  each  of 
the  radiation  references  as  well  as  the  terrain  (see  "ground  scan"), 
only  the  grayboay  references  apply  to  every  wavelength  band.  These 
graybody  references  (#1  and  #2  and  thermal  ambient)  serve  as  temperature 
calibration  sources  for  the  thermal  detectors  and  also  as  a  dark  level 
source  for  the  shorter-wavelength  non-thermal  detectors.  The  remaining 
sources  (lamp  and  sky)  serve  as  references  for  the  non-thermal  bands 
(as  shown).  For  indexing  purposes,  synchronization  references  are 
generated  by  the  scanner  and  recorded  with  the  video  signals.  The 
marker  pulse  refers  to  the  scan  position  relative  to  the  internally 
mounted  radiation  references;  the  roll-stabilized  pulse  refers  to 
ground  scan  nadir  with  aircraft  roll  motion  removed. 

Table  2.1  lists  significant  parameters  of  the  M-7  scanner  system. 
The  scanner  views  the  terrain  during  90°  of  its  scan,  providing  an 
external  field  of  view  (FOV)  +  45°  from  nadir.  A  nominal  0.1°C  NEAT 
and  a  1%  NEAp  are  achieved.*  The  system  operates  at  either  of  two 
constant  scan  speeds  -  60  or  100  scans  per  second.  Electronic  band¬ 
width  is  tape-recorder-limited  to  a  range  of  dc  to  90  kHz.  Table  2.2 
identifies  those  detector  assemblies  currently  in  use  with  the  system. 
Where  there  is  a  choice  of  detectors,  the  first-listed  unit  is  the  one 
commonly  used. 

The  M-5  scanner,  shown  schematically  in  Figure  2.5,  is  a  double- 
ended  scanner  using  a  double  "axe-blade"  scanning  mirror  to  direct 
radiation  to  the  two  ends.  In  most  respects  it  is  the  same  as  the 
M-7  except  that  the  data  collected  from  the  two  ends  of  the  scanner 

*NEAt  =  Noise  Equivalent  change  in  Temperature 
NEAp  =  Noise  Equivalent  change  in  reflection 
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FIGURE  2.  A  SCANNER  VOLTAGE  OUTPUT  VERSUS  TIME 


TABLE  2.1 


M7  SCANNER  PERFORMANCE  CHARACTERISTICS 

12  Spectral  Bands  In  UVf  Visible  and  IR  Regions 
90°  External  FOV  (j^5°  from  nadir) 

2  mrad  Maximum  Spatial  Resolution 
0.1°C  Nominal  Thermal  Resolution 
1 %  Nominal  Reflectance  Resolution 
Five  Radiation  Reference  Ports 
5-inch  Diameter  Collector  Optics 
Scan  Rate  of  60  or  100  scans /sec 
DC  to  90  kHz  Electronic  Bandwidth 
Roll-Stabilized  Imagery 


TABLE  2.2 


Top  of  Aircraft  Fuselagi 


Fiber-Optic  Bundle* 


Photomultiplier*  (12)- 


'  '  ; 

.  Fiber-Optic  Array 
•  Peimaging  Lent 


12-Channel 

Spectrometer 


Dispersing* 

/  Prisrir  ^ 


Entrance 
/Slit  ^ 


Collimating 
Lens  4 


Paraboloidal  Mkror 


Scanner  2;  End  A  "  " 


are  90°  out  of  phase.  Two  M-5  scanners  in  the  aircraft  were  used  to 
collect  data.  Essentially  the  same  detector  assemblies  were  used  with 
the  M-5  as  are  used  with  the  M-7  system.  One  M-5  scanner  was  used  to 
collect  multlspectral  data  in  the  visible  and  near  infrared.  Calibration 
lamps  were  added  in  the  scanner  housing  which  did  not  restrict  the 
f ield-of-view  below  the  aircraft.  Another  M-5  scanner  was  used  to  col¬ 
lect  thermal  data,  and  "hot"  and  "cold"  graybody  thermal  references 
were  added  at  the  scanner  aperture  so  the  thermal  IR  channels  could  be 
calibrated.  These  thermal  plates  extended  into  the  f ield-of-view  of 
the  scanner  below  the  aircraft  so  that  the  video  was  limited  to  approx¬ 
imately  +  20°  from  nadir. 
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PROCEDURES  AND  SELECTED  EXAMPLES 

The  multispectral  scanner  data  selected  for  terrain  backgrounds 
statistical  analysis  on  this  program  were  recorded  on  analog  magnetic 
tape,  and  these  analog  tapes  were  converted  to  high  density  digital  tapes. 
The  general  pre-processing  procedure  used  to  create  calibrated  data  tapes 
is  discussed  in  Section  3.1  with  a  detailed  discussion  of  the  procedures 
used  for  each  scene  being  presented  in  Appendix  I.  The  statistical 
measures  derived  from  the  data  are  described  in  Section  3.2;  point 
statistics  in  Section  3.2.1  and  correlative  statistics  in  Section  3.2.2. 
Three  types  of  correlative  statistics  are  discussed.  These  include 
spectral  correlations  between  pairs  of  bands;  area/ intensity  statistics 
and  the  development  of  the  equivalent  elliptical  area  concept;  and  one- 
and  two-dimensional  Wiener  (power)  spectra.  In  Section  3.3  an  example 
is  shown  in  which  the  original  background  scene  is  actually  replaced 
by  the  equivalent  elliptical  areas  and  most  of  the  spatial  information 
is  preserved. 

Examples  of  each  type  of  statistical  data  are  presented  in  Section 
3,  with  the  complete  summary  of  terrain  backgrounds  statistical  data 
included  as  Appendix  II. 

3.1  PREPROCESSING  OF  SCANNER  DATA 

Output  from  the  ERIM  scanners  was  converted  to  high-density- 
digital  tapes  in  which  each  data  value  is  recorded  as  an  8-bit  integer 
ranging  in  value  from  0  to  255.  Each  channel  (wavelength  band)  is 
recorded  on  a  separate  tape  channel  and  the  amplifier  gains  adjusted 
so  that  all  data  values  fall  in  the  0  -  255  range.  On  this  tape  the 
scan  lines  for  a  single  M-7  scanner  channel  or  reflective  IR  M-5  scanner 
channel  consists  of  790  data  points  of  which  the  first  646  are  scene 
elements  (pixels)  covering  the  range  -45°  to  +45°  with  respect  to  nadir 
while  the  remaining  144  points  are  calibration  values,  24  data  points 
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for  each  of  6  calculation  sources.  For  the  thermal  M-5  scanner  channels 
the  scanner  elements  cover  pixels  162  to  484  with  pixels  1  -  161  and 
485  -  646  being  used  for  the  graybody  thermal  sources. 

Before  the  scanner  data  may  be  used  to  generate  scene  statistics 
some  pre-processing  is  required: 

1)  The  high  density  digital  tapes  must  be  converted  to 
computer  compatible  tapes: 

2)  Each  channel  must  be  calibrated  in  temperature  or 
radiance  using  the  calibration  data  at  the  end  of  each 
scan  line; 

3)  Averaging  was  employed  to  reduce  oversampling  and  to 
equalize  any  differences  in  the  f ields-of-view  of  the 
various  detectors;  and 

4)  A  set  of  calibrated,  formatted  tapes  must  be  generated. 

The  first  of  these  processes,  conversion  of  high  density  tapes  to 

computer  compatible  tapes,  has  been  accomplished  using  the  conversion 
facilities  at  Bendix  Aerospace  Systems  Division  in  Ann  Arbor  and  the 
MIDAS  computer  system  at  ERIM.  The  result  of  this  conversion  process 
was  a  low  density  (800  BPI)  tape  for  each  of  the  scenes  desired.  Using 
these  tapes,  the  remainder  of  the  pre-processing  (2,3,  and  4)  was 
performed  using  a  computer  code  written  for  the  University  of  Michigan's 
AMDAHL  computer  system.  It  is  the  set  of  calibrated  tapes  generated 
from  this  code  that  are  used  for  all  image  processing. 

The  data  values  appearing  on  the  calibrated  tapes  are  themselves 

integers,  ranging  from  0  to  255,  but  these  integers  have  been  modified 

so  that  a  linear  relationship  exists  between  them  and  the  apparent 

scene  radiance  or  temperature.  The  data  in  the  near  IR  channels  was 

2 

converted  to  equivalent  radiance  in  uw/cm  •sr*ym).  The  equivalent 
radiance  is  the  value  of  the  spectral  radiance  at -the  center  wave¬ 
length  of  the  filter  produced  by  a  2850  K  NBS  lamp  source  filling  the 
sensor  aperture  and  giving  the  same  detector  response.  The  data  in- 
the  thermal  IR  channels  was  converted  to  apparent  temperature  in  degrees 
Kelvin.  The  apparent  temperature  is  the  temperature  of  a  blackbody 
filling  the  sensor  aperture  producing  the  same  detector  response. 


'  -PMWii-  u  i  Mwiwipuwui  m  jii 


A  table  of  apparent  temperatures  and  their  corresponding  band 
radiances  is  included  as  Table  1-1  in  Appendix  I.  The  radiance  (or 
temperature)  is  recovered  from  the  integer  data  values  using  multi¬ 
plicative  and  additive  factors  recorded  in  the  tape  header  record  for 
each  channel.  These  "mult"  and  "add"  factors  were  determined  from  the 
calibration  sources  appearing  at  the  end  of  each  scan  line  assuming 
a  linear  relationship  between  radiance  and  detector  output*.  For 
channels  calibrated  in  radiance,  the  integer  values  on  the  calibrated 
tape  differ  from  those  on  the  original  computer  compatible  tapes  by  at 
most  a  zero  level  correction.  For  channels  calibrated  in  temperature, 
however,  the  new  integers  are  distinctly  different  from  the  original 
values  since  those  on  the  calibrated  tape  are  linear  in  temperature 
while  those  on  the  original  tape  are  proportional  to  radiance. 

Averaging  was  employed  to  reduce  oversampling  and  to  equalize  any 
differences  in  the  f ields-of-view  of  the  various  detectors.  The  degree 
of  oversampling  was  generally  largest  in  the  along-track  direction 
because  of  the  constant  (60/sec)  scan  rate.  The  degree  of  over  sampling 
in  the  cross-track  direction  was  small  in  almost  all  of  the  data  because 
a  3.8  x  lO'Vsec.  sampling  rate  was  used  to  produce  a  sample  every  2.5 
mr  to  correspond  to  the  2.5  mr  cross-track  f ield-of-view  of  most  of  the 
detectors.  The  technique  that  has  been  used  for  averaging  lines  to 
reduce  oversampling,  rather  than  dropping  an  appropriate  number  of 
alternate  lines,  was  developed  for  two  reasons:  1)  to  improve  the 
signal-to-noise  ratio  in  the  resulting  data,  and  2)  to  equalize  any 
differences  in  the  along-track  fields-of-view  of  the  various  detectors 
for  calculating  the  spectral  correlation  coefficients.  The  same  averaging 
technique  can  be  used  to  reduce  any  overlap  in  the  cross-track  direction 
and  to  equalize  any  differences  in  the  cross-track  fields-of-view  of 
the  various  detectors,  but  cross-track  averaging  was  not  found  to  be 
necessary  in  any  of  the  data.  Figure  3-1  is  a  schematic  representation 
of  the  procedure  used*.  In  this  example,  D  is  the  largest  along-track 
f ield-of-view  of  the  channels  and  d  the  f ield-of-view  of  the  channel 

*  See  Appendix  I  for  a  detailed  description  of  the  procedures  used. 
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INPUT  OUTPUT 


FIGURE  3-1  SCHEMATIC  REPRESENTATION  OF  THE  PROCEDURES  USED  FOR  FIELD-OF-VIEW  AVERAGING 


being  averaged.  The  averaging  was  done  by  summing  the  radiances  of 
the  scan  line  times  their  overlap  with  the  desired  output  line  and 
dividing  by  the  sum  of  the  overlap  factors.  If  the  scan  lines  for  the 
largest  field-of-view  were  themselves  overlapped,  non-overlapped 
f ields-of-view  were  generated  by  taking  as  output  lines  those  with 
f ields-of-view  D  each  of  which  was  displaced  by  D  as  in  the  right 
side  of  Figure  3-1.  Overlap  factors  were  then  determined  in  an 
identical  manner  between  the  original  scan  lines  with  field-of-view 
D  and  this  set  of  non-overlapped  output  lines. 

After  calibration  and  field-of-view  averaging,  a  new  data  tape 
was  generated.  To  be  compatible  with  existing  data  processing  systems, 
this  tape  was  written  in  ERIM-7094  format  which  consists  of  36-bit 
words  each  of  which  contain  4  data  values.  The  individual  scan  lines 
were  written  with  the  channels  interleaved  with  646  9-bit  data  points 
per  channel  per  scan  line.  The  "mult"  and  "add"  factors  required 
for  calibration  of  the  data  were  written  in  the  tape  header  record 
along  with  necessary  format  information.  These  tapes  could  then  be 
used  directly  as  input  to  new  or  existing  statistics  generation  programs. 

3.2  STATISTICS  DEFINITIONS 

Several  sets  of  statistics  have  been  generated  for  each  of  the 
chosen  scenes.  These  may  broken  down  into  two  groups: 

1)  Point  statistics:  Those  defined  by  individual  data 
points  in  a  single  channel;  and 

2)  Correlative  and  area  statistics:  Those  requiring  cal¬ 
culation  of  correlation  effects  for  a  scene  either  in 

a  single  channel  or  between  channels. 

3.2.1  POINT  STATISTICS:  MEANS  AND  STANDARD  DEVIATIONS 

The  point  statistics  generated  were  the  mean  and  standard  deviation 
for  each  channel  and  a  histogram  of  the  data  value  distributions  of 
these  channels.  To  determine  the  degree  of  homogenity  that  existed, 
the  total  scene  was  broken  down  into  sub-areas  and  the  point  statistics 
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generated  for  these  sub-areas  as  well  as  for  the  total  area.  In 

practice,  the  sub-area  statistics  were  generated  first  and  the  total 

area  statistics  derived  from  them. 

The  mean  value  for  sub-area  n  in  Channel  J  (x(J)  )  was  evaluated 

0 

using 

Nn 

Z  *(J>i  u> 

n  i«i 


where  x(J) .  is  the  data  value  of  pixel  i  in  Channel  J  and  N  is  the 

i  n 

total  number  of  data  points  in  the  sub-area.  Using  the  same  notation 
the  standard  deviation,  0(1)^,  is  given  by 


The  corresponding  total  scene  values  in  Channel  J  are  then  related  to 
those  of  the  sub-areas  through  the  following  relations 

x(J)  =  ^  Z  N  x(J) 

n  *-  n  n 
n 


and 


(3) 


o2(  J) 


(«4r) 


_ 2~ ]  _ 2 

+  -  x(J) 


where  N  is  the  number  of  points  in  the  total  scene.  Typical  examples 
of  these  statistics,  generated  from  the  Flint-1  data,  are  given  in 
Table  3-1.  In  this  scene  the  sub-areas  had  on  the  order  of  77  000 
data  points  while  the  total  scene  had  approximately  464,000  points. 

At  the  same  time  that  these  statistics  were  generated,  the  number 
of  data  points  having  each  of  the  possible  data  values  (0  to  255)  was 
tabulated  and  this  tabulation  used  to  generate  histograms  for  the  scene. 
Representative  sub-areas  and  total  area  histograms  for  the  9.3  to  11.7  ym 
channel  of  Flint-1  are  shown  in  Figures  3-2  and  3-3  respectively. 
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3.2.2  CORRELATIVE  AND  AREA  STATISTICS 

Three  sets  of  correlative  statistics  were  generated  for  the  scenes: 
wavelength  correlations,  area  statistics,  and  spatial  correlations  or 
Weiner  spectra.  Because  of  the  common  factors  involved,  the  wavelength 
correlations  were  evaluated  at  the  same  time  as  the  means  and  standard 
deviations  while  the  area  statistics  and  Weiner  spectra  required  separate 
computer  codes. 

WAVELENGTH  CORRELATIONS 

The  wavelength  correlations  were  determined  for  all  pairs  of 
channels  except  for  those  collected  at  opposite  ends  of  the  M-5  which 
were  90°  out  of  phase.  The  data  for  the  1.0-1. 4,  1.5-1. 8,  and  2. 0-2. 6  ym 
channels  were  collected  using  a  three  element  InAs  detector,  the  seg¬ 
ments  of  which  look  at  slightly  different  ground  areas  at  any  given  time. 
The  1.0-1. 4,  2. 0-2. 6,  and  4. 5-5. 5  ym  data  were  collected  with  a  three- 
element  InSb  detector.  The  three  channel  data  were  brought  into  regis¬ 
tration  at  the  nadir  by  applying  appropriate  time  delays  to  the  leading 
channels  (time  delays  long  enough  to  correct  for  the  90°  phase  difference 
between  the  two  ends  of  the  M-5  scanner  could  not  be  generated) . 

Complete  registration  could  not  be  maintained  across  the  entrie  scan 
line  because  as  the  scan  angle  is  increased  the  projection  of  the  three 
detector  array  rotates  about  the  center  element  until  at  a  scan  angle 
of  90°  the  detector  array  is  aligned  in-track  rather  than  cross-track. 
However,  the  element-to-element  registration  is  50%  or  better  for  angles 
between  +  20°  of  nadir.  Since  wavelength  correlations  were  calculated 
at  the  same  time  as  the  means  and  standard  deviations,  the  sub-areas 
discussed  above  were  chosen  so  that  one  sub-area  always  covered  the 
+  20°  field.  Wavelength  correlations  were  generated  for  all  sub-areas 
but  because  of  the  lack  of  complete  registration  those  for  the  +  20° 
are  most  meaningful. 

The  wavelength  correlation  between  two  channels,  J  and  K,  is  defined 
in  terms  of  a  correlation  coefficient  which  is  the  ratio  of  the  covar¬ 
iance  of  the  two  channels  to  the  product  of  their  standard  deviations: 
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C0R(J  K)  =,  COV^JjK) 

cor(j,k;  CT(J)a(K)  (4) 

where  COR(J,K)  is  Che  correlation  coefficient  for  the  two  channels, 

a(J)  and  a(K)  their  standard  deviations,  and  COV(J,K),  the  J-K  element 

/ 

of  the  covariance  matrix  for  a  given  sub-area  ,  is  defined  as 

Nn 

C0V<J-K>n  L  (5) 

where  all  symbols  are  the  same  as  those  of  Equations  1  and  2.  The 

reason  for  the  simultaneous  evaluation  of  Equations  2  and  5  is  obvious: 

2 

the  square  of  the  standard  deviation  a  (J)  is  the  autocovariance  or  the 
covariance  of  a  given  channel  with  itself  COV(J,J). 

An  example  of  the  correlation  coefficients  for  the  infrared  chan¬ 
nels  of  Flint-1  are  given  in  Table  3-2. 

AREA/ INTENSITY  STATISTICS 

Since  means,  standard  deviations,  and  histograms  do  not  give  any 
information  about  the  positions  of  data  values  in  the  scene  or  possible 
clustering  of  these  values,  area  statistics  have  been  generated  by  deter¬ 
mining  contiguous  regions  of  the  scene  for  which  the  enclosed  points  had 
values  greater  than  some  prescribed  threshold  [Reference  3].  Once 
these  regions  had  been  defined,  their  geometric  centroids,  areas,  and 
second  moments  were  determined  and  these  parameters  used  to  define  a 
set  of  equivalent  elliptical  areas.  The  output  from  this  procedure  was, 
for  each  threshold  level,  a  tabulation  of  centroids,  tilt  angles  and 
semi -major  and  minor  axes  defining  the  equivalent  ellipses. 

The  semi -major  and  semi-minor  axes  of  the  equivalent  ellipses  were 

taken  coincident  with  the  principal  axes  of  the  regions  they  represented, 
c  c 

If  I  and  I  are  the  second  moments  of  a  region  about  the  centroid  and 
c  x  y 

I  its  product  of  inertia,  the  angle  of  rotation  of  the  principal  axes 
xy 

relative  to  a  fixed  x-y  coordinate  system  is  given  by 

[3]  W.  G.  Burge,  W.  L.  Brown,  A  study  of  Waterfowl  Habitat  in  North 

Dakota  Using  Remote  Sensing  Techniques,  Report  No.  2771-7-F,  Willow 
Run  Laboratories,  University  of  Michigan,  Ann  Arbor,  July  1970. 
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TABLE  3-2 

CORRELATION  COEFFICIENTS  FOR  THE  TOTAL  AREA  OF  FLINT-1 


Spectral 


Correlation 

1.0-1. 4  um 

1.5-1. 8  um 

2. 0-2. 6  um  9.3-11.7  um 

1.0-  1.4 

pm 

1.000 

1.5-  1.8 

pm 

0.392 

1.000 

2.0-  2.6 

pm 

0.303 

0.603 

1.000 

9.3-11.7 

pm 

-0.455 

0.048 

0.177  1.000 
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where  a  is  in  radians.  The  second  moments  of  the  region  about  the 
principal  axes  are  then  given  by 


cos  2a  +  I  sin  2a 
xy 


(7) 


Since  the  area  and  the  second  moments  of  an  arbitrary  region  can¬ 
not  be  simultaneously  matched  using  an  ellipse,  equality  was  demanded 
between  the  geometric  areas  and  the  ratios  of  the  principal  moments.  If 
the  semi-major  and  semi-minor  axes  of  the  ellipse  are  defined  as  a  and  b 
respectively  and  the  a-axis  is  aligned  with  the  larger  of  the  principal 
axes  of  the  region,  these  equalities  give 


(8) 


where  A  is  the  area  of  the  region  and  1^  and  I 2  are  the  moments  given 
by  Equation  7  with  ^  >  I2<  If  Iyt  of  Equation  7  is  the  larger  of 
the  moments,  the  a-axis  of  the  ellipse  is  oriented  at  an  angle  a  with 
respect  to  the  fixed  x-axis;  if  I  ,  is  larger,  the  a-axis  is  at  an  angle 
a  +  ir/2.  An  example  of  the  equivalent  elliptical  areas  generated  for  the 
9.3  -  11.7  pm  channel  of  Flint-1  is  shown  in  Figures  3-4  through  3-7. 

Each  of  these  figures  represents  the  areas  found  for  data  values 
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FIGURE  3-4  EQUIVALENT  ELLIPTICAL  AREAS 

FOR  THE  9.3  -  11.7  ym  CHANNEL 
OF  FLINT-1  WITH  A  1.5o  THRESHOLD 
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RE  3-6  EQUIVALENT  ELLIPTICAL  AREAS  F. 
THE  9.3  -  11.7  vm  CHANNEL  OF 
FT.TNT-1  WTTH  A  2.5cr  THRESHOLD 


exceeding  the  specified  threshold  (1.5o,  2.0a,  2.5a,  3.0 a)  where  a  is 
the  standard  deviation  (3.1°K)  and  the  threshold  is  defined  as  the  mean 
(294°K)  plus  the  specified  multiple  of  a. 

In  addition  to  generating  the  parameters  required  for  definition 
of  the  equivalent  ellipses,  the  perimeters  of  these  regions  have  been 
tabulated  and  shape  factors  evaluated.  These  shape  factors  (SF)  were 
defined  as 


SF  -  — - —  (9) 

2/iT  fk 

where  P  is  the  perimeter  of  the  region,  A  its  area,  and  the  factor  1/2/ir 

included  so  that  a  circular  area  gives  a  shape  factor  of  unity. 

For  all  regions  found,  tabular  histograms  have  been  formed  sorting 

the  regions  by  area,  perimeter,  and  shape  factor.  Examples  of  these 

histograms,  again  for  the  9.3-11.7  um  channel  of  Flint-1,  are  given  in 

Table  3-3  through  3-5  corresponding  to  a  1.0a  threshold.  Tabulations  with 

increments  in  area  and  perimeter  down  to  the  resolution  of  the  scanner 

2 

data,  approximately  1  m  ,  are  available. 

WIENER  (POWER)  SPECTRA 

The  last  of  the  correlative  statistics  generated  were  Wiener 
spectra  which  gave  information  about  the  spatial  frequency  content  of 
the  images.  The  Weiner  spectrum  for  a  stationary  process  is  defined 
as  the  Fourier  transform  of  the  autocorrelation  function.  In  one 
dimension  this  is 


S(k  )  =  f  exp  (-2rcik  x)R(x)dx  (10) 

x  -4  x 

where  S(k^)  is  the  Wiener  or  power  spectrum,  k^  the  spatial  frequency 
and  R(x)  the  autocorrelation  function.  R(x)  is  itself  defined  as 
the  expectation  value  of  the  product  of  scene  data  values  times  the 
corresponding  values  for  the  scene  when  displaced  by  x: 

R(x)  -  E  (f(X)f(X  +  x)}  (11) 
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TABLE  3-3 

AREA  DISTRIBUTION  FOR  RECOGNIZED  REGIONS 
IN  THE  9.3-11.7  urn  CHANNEL  OF  FLINT-1  WITH  A  1 o  THRESHOLD 


SQUARE 

METERS 

FREQUENCY 

0.0 

TO 

100.0 

2047 

100.0 

TO 

200.0 

78 

200.0 

TO 

500.0 

26 

500.0 

TO 

1000.0 

8 

1000.0 

TO 

1500.0 

2 

1500.0 

TO 

2000.0 

0 

2000.0 

TO 

2500.0 

0 

2500.0 

TO 

3000.0 

0 

3000.0 

TO 

4000.0 

1 

4000.0 

TO 

5000.0 

0 

5000.0 

TO 

6000.0 

0 

6000.0 

TO 

8000.0 

1 

8000.0 

TO 

10000.0 

0 

10000.0 

TO 

15000.0 

0 

15000.0 

TO 

20000.0 

0 

20000.0 

TO 

40000.0 

0 

40000.0 

TO 

80000.0 

0 

80000.0 

TO 

160000.0 

0 

OVER 

160000.0 

0 

TABLE  3-4 

PERIMETER  DISTRIBUTIONS  FOR  RECOGNIZED  REGIONS 
IN  THE  9.3-11.7  ura  CHANNEL  OF  FLINT-1  WITH  A  1  a  THRESHOLD 


METERS 

FEET 

FREQUENCY 

0 

TO 

50 

0 

TO 

164 

1991 

50 

TO 

100 

164 

TO 

328 

128 

100 

TO 

150 

328 

TO 

492 

16 

150 

TO 

200 

492 

TO 

656 

15 

200 

TO 

250 

656 

TO 

820 

5 

250 

TO 

300 

820 

TO 

984 

3 

300 

TO 

350 

984 

TO 

1148 

1 

350 

TO 

400 

1148 

TO 

1312 

1 

400 

TO 

500 

1312 

TO 

1640 

1 

500 

TO 

600 

1640 

TO 

1968 

0 

600 

TO 

700 

1968 

TO 

2296 

0 

700 

TO 

800 

2296 

TO 

2624 

0 

800 

TO 

900 

2624 

TO 

2952 

1 

900 

TO 

1000 

2952 

TO 

3280 

1 

1000 

TO 

1200 

3280 

TO 

3937 

0 

1200 

TO 

1400 

3937 

TO 

4593 

0 

1400 

TO 

1600 

4593 

TO 

5249 

0 

1600 

TO 

2000 

5249 

TO 

6561 

0 

OVER 

2000 

OVER 

6561 

0 

42 


H 


/ 


TABLE  3-5 

SHAPE  FACTOR  DISTRIBUTIONS  FOR  RECOGNIZED  REGIONS 
IN  THE  9.3-11.7  um  CHANNEL  OF  FLINT-1  WITH  A  1 o  THRESHOLD 


SHAPE  FACTOR  FREQUENCY 


0.0  TO 

1.0 

11 

1.0  TO 

1.1 

2 

1.1  TO 

1.2 

447 

1.2  TO 

1.3 

303 

1.3  TO 

1.4 

513 

1.4  TO 

1.5 

280 

1.5  TO 

1.6 

227 

1.6  TO 

1.7 

133 

1.7  TO 

1.8 

59 

1.8  TO 

1.9 

50 

1.9  TO 

2.0 

34 

2.0  TO 

2.2 

43 

2.2  TO 

2.4 

21 

2.4  TO 

2.6 

9 

2.6  TO 

2.8 

7 

2.8  TO 

3.0 

5 

3.0  TO 

3.5 

15 

3.5  TO 

4.0 

3 

OVER 

4.0 

1 

43 


where  E  represents  the  expectation  value  of  the  argument.  The  Wiener  spec¬ 
trum  may  be  evaluated  without  first  determining  the  autocorrelation 
function  if  the  integral 


00 

j  |xR(x)|< 


is  bounded,  which  is  usually  the  case  for  non-periodic  data  with  zero 
mean.  In  this  case  it  may  be  shown  that  [Reference  A] 


S(kx)  =  lim  (x2 


-  X;L)->»  J2  f (x)  exp(-2irikxx)dx 


so  that  the  Wiener  spectrum  is,  in  the  limit,  the  modulus  squared  of 
the  Fourier  transform  of  the  scene.  Written  in  terms  of  the  discrete 
Fourier  transform,  this  equation  becomes 

2 

Ax  NfX 

S(j)  -  ft - T)  f(£)  exp  <-2"ij*/N)  (1< 


where  N  is  the  total  number  of  points  being  transformed,  i  is  the  square 
root  of  -1,  and  the  spatial  frequency  kx>  evaluated  only  at  integer 
values  of  j ,  is  given  by 


\e  m  - 

x  NAx 


where  Ax  is  the  displacement  between  successive  data  points. 

It  may  be  seen  from  Equation  14  that  the  Wiener  spectrum  S(j) 
is  symmetric  about  j  ■  N/2  since 

exp  (-2vi(N  -  j)£/l 0“  exp  (+2irij£/N) 


A.  Papoulis,  Probability,  Random  Variables,  and  Stochastic  Processes, 
McGraw-Hill  Book  Company,  New  York,  1965. 


llgWWI 


I 


Hence 


S(N  -  j)  -  S(j)  for  j  -  1,  2 . N/2 


and  only  che  first  half  of  the  Wiener  spectrum  needs  to  be  evaluated. 
The  calculated  frequency  range  is  then 


JL  <  k  <  _A 

NAx  -  x  -  2  Ax 


The  two  dimensional  Wiener  spectrum  is  defined  in  an  analgous  manner 
except  that  two  transforms  are  required,  one  for  each  dimension: 


S(j’  “  (bTt)  (iT^l)  Z!  exp  <*2irijn/N)* 


M-l  2 

exp  (-2iri-£m/M)f  (n,m) 

m*0 


where  the  spatial  frequencies  in  the  two  dimensions  are  given  in  terms 
of  the  integers  j  and  L  by 


k*  NAx 


One  dimensional  Wiener  spectra  were  evaluated  both  cross-track  (along 
scan  line)  and  in-track  (parallel  to  the  aircraft  flight  paths) .  These 
are  average  spectra  since  Equation  14  was  used  to  transform  individual 
"lines"  of  data  and  the  transforms  have  been  averaged  over  a  number  of 
"lines".  Since  the  Fourier  transform  algorithms  were  substantially  less 
expensive  if  the  number  of  points  transformed  was  a  power  of  2,  only  512 
points  were  used  for  both  the  cross-track  and  in-track  "lines"  of  data. 
The  middle  512  scan  lines  were  averaged  to  obtain  the  cross-track  Wiener 
spectra;  25  along-track  lines  equally  spaced  across  the  image  were  aver¬ 
aged  to  obtain  the  in-track  Wiener  spectra.  Examples  of  in-track  and 
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cross-track  Wiener  spectra  for  the  9.3-11.7  pm  channel  of  Flint-1  are 

2 

shown  in  Figures  3-8  and  3-9.  The  units  of  power  density  are  (K)  /cycle 

2  2 

per  meter  for  the  9.3-11.7  pm  band  and  (pw/cm  *sr*pm)  /cycle  per  meter 

for  the  1.0-1. 4,  1.5-1. 8  and  2.0-2. 6  pm  bands. 

Representative  two-dimensional  Wiener  spectra  were  evaluated  only 

for  the  Flint-1  scene.  Figure  3-10  shows  a  plot  of  the  two-dimensional 

spectrum  for  the  9.3-11.7  pm  channel.  The  units  of  power  density  for 

2  2 

the  two-dimensional  Wiener  spectrum  are  (K)  / (cycle  per  meter)  .  As 
for  the  one-dimensional  spectra,  the  two-dimensional  spectra  were 
evaluated  using  512  data  points  in  each  of  the  dimensions.  The  utility 
of  the  two-dimensional  Wiener  spectra  to  the  system  designer  has  not 
been  firmly  established  so  that  only  limited  two-dimensional  Wiener 
spectra  processing  has  been  performed  to  date. 


FIGURE  3-10  TWO-DIMENSIONAL  WIENER  SPECTRUM  FOR  THE  9.3  -  11.7  Pm  CHANNEL  OF  FLINT- 


3.3  MODELING  THE  BACKGROUND  SCENE 

Although  the  statistics  discussed  above  give  information  about  the 
particular  scene,  they  do  not  provide  a  means  of  modeling  the  scene  or 
of  generating  "typical"  scenes.  A  preliminary  modeling  procedure  has 
been  developed  using  the  concept  of  elliptical  areas  as  discussed  in 
connection  with  the  area  statistics  above.  Basically  the  scene  is  bro¬ 
ken  down  into  elliptical  areas  at  several  threshold  levels,  as  for  the 
area  statistics,  and  each  pixel  within  an  area  is  given  a  value  equal  to 
the  threshold.  When  this  is  done  for  a  number  of  thresholds,  the  result 
is  a  set  of  three-dimensional  elliptical  cylinders  whose  area  is  the 
geometric  area  of  the  original  scene  element  and  whose  height  corresponds 
to  that  of  the  threshold.  These  cylinders  when  superimposed  give  a 
pseudo-image  whose  spatial  characteristics  approximate  those  of  the 
original  image  for  many  sensor  systems  performance  analyses.  The  advan¬ 
tage  of  this  scene  is  that  the  individual  areas  may  be  moved  at  random 
throughout  the  scene  and  "typical"  scenes  may  be  generated  by  defining 
distributions  for  the  various  sets  of  elliptical  volumes.  Parameters 
defining  a  typical  scene  could  then  be  distributions  of  geometric  areas 
for  given  thresholds  and  a  joint  density  distribution  giving  the  prob¬ 
ability  that  if  a  cylinder  of  a  given  size  and  threshold  is  present, 
a  smaller  cylinder  of  a  higher  threshold  and  size  will  be  coincident 
with  it. 

At  the  present  time,  one  pseudo-image  has  been  generated  for  the 

9.3  -  11.7  pm  channel  of  Flint-1.  Although  area  distributions  are 
generated  at  the  time  the  various  areas  are  found,  no  attempt  has  been 
made  to  determine  the  joint  probabilities  *  Figures  3-11  and  3-12  show 
artificial  color  representations  of  the  original  image  and  the  psuedo- 
image  respectively.  For  this  specific  case  the  threshold  levels  have 
been  taken  arbitrarily  at  la,  2a,  3a,  and  4a  (a  •  3.1  K)  above  the 
mean  (294  K)  and  these  levels  represented  ty  the  colors  blue,  green, 
yellow,  and  red  respectively.  The  area  distributions  are  shown  in 
Tables  3-6  through  3-9. 
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TABLE  3-6 

AREA  DISTRIBUTIONS  FOR  RECOGNIZED  REGIONS  IN  THE 
REAL  AND  PSEUDO-IMAGES  OF  THE  9.3-11.7  gm  CHANNEL  OF  FLINT-1. 

THRESHOLD  LEVEL  =  1  a 


SQUARE 

METERS 

FREQUENCY 

0.0 

TO 

100.0 

2047 

100.0 

TO 

200.0 

78 

200.0 

TO 

500.0 

26 

500.0 

TO 

1000.0 

8 

1000.0 

TO 

1500.0 

2 

1500.0 

TO 

2000.0 

0 

2000.0 

TO 

2500.0 

0 

2500.0 

TO 

3000.0 

0 

3000.0 

TO 

4000.0 

1 

4000.0 

TO 

5000.0 

0 

5000.0 

TO 

6000.0 

0 

6000.0 

TO 

8000.0 

1 

8000.0 

TO 

10000.0 

0 

10000.0 

TO 

15000.0 

0 

15000.0 

TO 

20000.0 

0 

20000.0 

TO 

40000.0 

0 

40000.0 

TO 

80000.0 

0 

80000.0 

TO 

160000.0 

0 

OVER 

160000.0 

0 

TOTAL  NUMBER  OF  REGIONS  =  2163 
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TABLE  3-7 


AREA  DISTRIBUTIONS  FOR  RECOGNIZED  REGIONS  IN  THE 
REAL  AND  PSEUDO-IMAGES  OF  THE  9.3-11.7  um  CHANNELS  OF  FLINT-1. 

THRESHOLD  LEVEL  =2  a. 


SQUARE 

METERS 

FREQUENCY 

0.0 

TO 

100.0 

1110 

100.0 

TO 

200.0 

10 

200.0 

TO 

500.0 

4 

500.0 

TO 

1000.0 

2 

1000.0 

TO 

1500.0 

0 

1500.0 

TO 

2000.0 

0 

2000.0 

TO 

2500.0 

0 

2500.0 

TO 

3000.0 

0 

3000.0 

TO 

4000.0 

0 

4000.0 

TO 

5000.0 

0 

5000.0 

TO 

6000.0 

1 

6000,0 

TO 

8000.0 

0 

8000.0 

TO 

10000.0 

0 

10000.0 

TO 

15000.0 

0 

15000.0 

TO 

20000.0 

0 

20000.0 

TO 

40000.0 

0 

40000.0 

TO 

80000.0 

0 

80000.0 

TO 

160000.0 

0 

OVER 

160000.0 

0 

TOTAL  NUMBER  OF  REGIONS  -  1127 


54 


TABLE  3-8 


AREA  DISTRIBUTIONS  FOR  RECOGNIZED  REGIONS  IN  THE 
RFAL  AND  PSEUDO-IMAGES  OF  THE  9.3-11.7  um  CHANNEL  OF  FLINT-1. 

THRESHOLD  LEVEL  =3  a. 


SQUARE 

METERS 

FREQUENCY 

0.0 

TO 

100.0 

97 

100.0 

TO 

200.0 

1 

200.0 

TO 

500.0 

2 

500.0 

TO 

1000.0 

1 

1000.0 

TO 

1500.0 

0 

1500.0 

TO 

2000.0 

0 

2000.0 

TO 

2500.0 

0 

2500.0 

TO 

3000.0 

0 

3000.0 

TO 

4000.0 

1 

4000.0 

TO 

5000.0 

0 

5000.0 

TO 

6000.0 

0 

6000.0 

TO 

8000.0 

0 

8000.0 

TO 

10000.0 

0 

10000.0 

TO 

15000.0 

0 

15000.0 

TO 

20000.0 

0 

20000.0 

TO 

40000.0 

0 

40000.0 

TO 

80000.0 

0 

80000.0 

TO 

160000.0 

0 

OVER 

160000.0 

0 

TOTAL  NUMBER  OF  REGIONS  -  102 


TABLE  3-9 

AREA  DISTRIBUTIONS  FOR  RECOGNIZED  REGIONS  IN  THE 
REAL  AND  PSEUDO-IMAGES  OF  THE  9.3-11.7  ym  CHANNEL  OF  FLINT-1. 

THRESHOLD  LEVEL  =  4  a. 


SQUARE 

METERS 

FREQUENCY 

0.0 

TO 

100.0 

6 

100.0 

TO 

200.0 

1 

200.0 

TO 

500.0 

0 

500.0 

TO 

1000.0 

1 

1000.0 

TO 

1500.0 

0 

1500.0 

TO 

2000.0 

1 

2000.0 

TO 

2500.0 

0 

2500.0 

TO 

3000.0 

0 

3000.0 

TO 

4000.0 

0 

4000.0 

TO 

5000.0 

0 

5000.0 

TO 

6000.0 

0 

6000.0 

TO 

8000.0 

0 

8000.0 

TO 

10000.0 

0 

10000.0 

TO 

15000.0 

0 

15000.0 

TO 

20000.0 

0 

20000.0 

TO 

40000.0 

0 

40000.0 

TO 

80000.0 

0 

80000.0 

TO 

160000.0 

0 

OVER 

160000.0 

0 

TOTAL  NUMBER  OF  REGIONS  -  9 
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DISCUSSION  OF  RESULTS 

The  complete  summary  of  statistical  measures  for  the  chosen  back¬ 
ground  scenes  are  included  as  Appendix  II  of  this  report.  Means,  standard 
deviations,  and  spectral  correlation  coefficients  are  tabulated,  and 
histograms  for  sub-areas  of  the  scene  as  well  as  for  the  total  scene 
are  presented. 

Some  of  the  more  obvious  features  follow  expected  patterns.  For 
example,  the  mean  values  of  radiances  in  the  reflective  near  IR  channels 
vary  widely,  due  in  part  to  the  different  background  reflectivities  in  the 
different  areas  and  in  part  to  variation  in  solar  elevation  angles  at 
the  time  the  data  were  collected  (giving  some  large  shadowed  areas) .  The 
Mill  Creek  data  were  collected  early  in  the  morning  and  the  radiance  levels 
are  low  and  shadows  clearly  evident.  The  standard  deviations  are  propor¬ 
tionately  larger  for  this  scene  because  of  these  shadows.  In  the  thermal 
channel,  9.3-11.7  Pm,  the  standard  deviation  about  the  mean  temperature 
is  smaller  in  the  Mill  Creek  data  than  in  other  scenes  as  is  expected 
early  in  the  morning.  The  spectral  correlation  coefficients  between 
pairs  of  near  IR  channel  data  are  high  in  the  Mill  Creek  data  and 
low  in  the  Baltimore  data.  This  is  indicative  of  a  wide  variety  of 
materials  in  the  Baltimore  scene,  some  with  reflectivities  that  increase 
and  others  that  decrease  with  wavelength,  and  of  the  large  shadowed  areas 
and/or  the  predominance  of  materials  in  the  Mill  Creek  scene  that  have  the 
same  relative  spectral  reflectance  characteristics.  The  spectral  correl¬ 
ations  between  the  reflective  and  thermal  IR  channels  are  low  as  expected. 

A  summary  of  these  total  area  statistics  is  shown  in  Table  4-1. 

The  histogram  data  reveal  the  very  non-Gaussian  characteristics  of 
most  of  the  terrain  backgrounds  data.  Especially  notable  is  the  wide 
distribution  of  the  Mill  Creek  data  with  its  large  areas  of  shadow. 
Histograms  derived  from  the  multispectral  scanner  data  of  Mill  Creek 
should  vary  dramatically  during  the  diurnal  cycle  because  of  the  topography 
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TABLE  4-1.  TOTAL  AREA  DATA  SUMMARY 


TABLE  4-1.  TOTAL  AREA  DATA  SUMMARY  (Continued) 
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and  this  factor  must  be  accounted  for  when  using  these  backgrounds 
data  for  sensor  performance  estimates  under  conditions  other  than 
those  for  which  the  data  were  collected.  The  histograms  for  the 
various  sub-areas  of  the  Baltimore  and  Mill  Creek  scenes  show  the 
largest  inhomogenie ty  of  the  background  scenes  analyzed. 

The  one-dimensional  power  spectra  for  the  1.0-1. 4  and  9.3-11.7  pm 
data  of  Flint-1  and  Mill  Creek  do  not  show  any  particularly  unique  or 
distinguishing  characteristics,  nor  do  the  two-dimensional  power  spectra 
for  the  Flint-1  data.  However,  there  are  very  pronounced  differences  in 
the  spatial  characteristics  of  these  background  scenes  which  may  be 
important  to  the  performance  of  a  sensor  in  detecting  or  tracking  a 
target  against  these  backgrounds.  The  area/intensity  statistics,  that 
is  the  occurrence  of  equivalent  elliptical  areas  representing  contiguous 
regions  of  the  background  above  a  threshold,  for  Flint-1  and  Mill  Creek 
are  quite  distinctive.  Especially  noticeable  in  the  thermal  infrared 
are  the  warm  roofs  of  houses,  many  of  which  exceed  two  standard  devia¬ 
tions  above  the  mean. 

The  area/intensity  statistics  do  preserve  many  of  the  character¬ 
istics  of  background  scenes  that  may  cause  false  alarms  in  many  sensor 
systems.  The  equivalent  elliptical  areas  above  each  threshold  setting, 
lo,  2 a,  3o,  and  4a  above  the  mean,  can  in  most  instances  be  identified 
with  the  specific  spatial  features  they  represent  in  the  original 
image.  A  simulation  of  the  actual  scene  with  these  equivalent  ellipti¬ 
cal  areas  was  shown  in  Section  3.  The  utility  of  the  area/intensity 
statistics  will  depend  on  whether  or  not  the  simulated  scene  adequately 
represents  the  actual  scene  for  sensor  performance  estimates,  and  whether 
or  not  sensor  performances  depends  specifically  on  the  placement  of 
these  equivalent  ellipses  within  the  scene.  This  subject  will  be  the 
basis  for  future  terrain  background  statistical  analyses. 


SUMMARY  AND  RECOMMENDATIONS 


Several  terrain  backgrounds  were  selected  and  statistics  were 
derived  for  these  backgrounds  from  calibrated  multispectral  scanner  data. 
Conventional  statistical  parameters  -  means,  standard  deviations,  histo¬ 
grams,  spectral  correlations,  and  power  spectra  -  are  reported.  In 
addition  a  program  has  been  initiated  for  developing  statistical  para¬ 
meters  of  terrain  backgrounds  particularly  relevant  to  the  sensor  de¬ 
signer's  problem  of  detecting  and  tracking  targets  with  low  probabil¬ 
ity  of  false  alarm  or  false  track.  These  statistical  parameters  are  the 
equivalent  elliptical  area/intensity  statistics  that  have  been  reported. 

It  Is  intuitively  obvious  that  these  area/intensity  statistics 
are  approximations  of  what  sensor  system  designers  need  to  estimate 
systems  performance  characteristics.  It  is  now  necessary  1)  to  demon¬ 
strate  that  these  statistics  are  adequate,  i.e.  that  the  performance 
of  sensors  against  the  actual  background  and  against  the  simulated 
background  is  essentially  the  same  and  that  the  salient  spatial  features 
of  the  background  have  been  preserved,  and  2)  to  demonstrate  that  ter¬ 
rain  backgrounds  can  be  classified  in  some  way,  e.g.  residential,  indus¬ 
trial,  mountainous,  etc.,  so  that  the  statistics  can  be  inferred  without 
recourse  to  actual  measurement.  Future  backgrounds  efforts  should  be 
directed  to  demonstrating  the  utility  of  the  area/ intensity  statistics 
for  terrain  backgrounds  and  to  expanding  these  efforts  to  include  sky 
and  cloud  backgrounds. 
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APPENDIX  I 

DETAILS  OF  CALIBRATION  AND  FIELD-OF-VIEW 
EQUALIZATION  PROCEDURES 

CALIBRATION  OF  SCANNER  DATA 

As  discussed  in  Section  3.1,  output  from  the  ERIM  scanners  is  in 
the  form  of  digital  tapes  in  which  the  data  values  are  represented  by 
integers  ranging  in  value  from  0  to  255.  At  the  end  of  each  scan  line 
a  set  of  calibration  sources  is  scanned  and  the  integers  observed  for 
these  sources  are  used  to  calibrate  the  data  in  apparent  radiance  or 
temperature.  The  resulting  calibrated  tape  also  represents  the  data  as 
integers  but  these  integers  are  adjusted  so  that  they  are  linear  in 
radiance  or  temperature  and  a  set  of  multiplicative  and  additive  factors 
are  used  to  convert  these  integer  values  to  the  appropriate  units. 

The  calibration  sources  scanned  once  per  scan  line  are:  controllable 
temperature  "hot"  and  "cold"  plates;  an  ambient  temperature  plate; 
a  visible-near-IR  lamp;  an  ultraviolet  lamp;  and  a  sun  sensor.  The 
visible-reflective  IR  channels  (approximately  0.4  pm  to  3.0  Pm)  are 
generally  calibrated  in  radiance  using  the  visible-near-IR  lamp  as  a 
radiance  standard  and  the  "ambient"  plate  as  a  dark  level  reference. 

If  a  linear  relationship  is  assumed  between  radiance  and  detector  output, 
the  apparent  radiance  of  the  target  (the  radiance  observed  at  the  scanner 
aperture)  is  given  by 

where 

-  The  integer  value  observed  for  the  target 

VA  -  The  integer  value  observed  for  the  ambient  plate 

-  The  integer  value  observed  for  the  lamp 
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The  radiance  of  the  lamp  at  the  center  of  the  bandpass 
of  the  channel  to  be  calibrated 


The  ambient  radiance  at  the  center  of  the  bandpass  of 
the  channel  to  be  calibrated 


and  is  the  apparent  radiance  of  the  target  which  is  related  to  the 
actual  target  radiance  as 


lt  “  ltvp  +  S 


(1-2) 


where  Tp  is  transmission  of  the  intervening  path  and  Lp  the  radiance 
of  this  path. 

The  lamp  radiance  in  Equation  1-1  is  obtained  from  a  calibration 
of  the  visible-near-IR  lamp  using  an  NBS  standard  lamp  while  the  ambient 
radiance  LA  is  taken  as  the  radiance  of  the  ambient  temperature  plate 
with  emissivity  t  given  by 


T  T  8b  .  / 1  .  -  BB 

LA  "  eALA  +  (1  "  eA)LA 


(1-3) 


where  the  first  term  is  the  plate  emittance  and  the  second  term  the 
surrounding  radiance  reflected  from  this  plate.  Since  the  ambient 
temperature  plate  and  its  surroundings  are  at  the  same  temperature, 
this  equation  simply  gives  L  equal  to  L“  or  the  radiance  of  a  black- 

CL  CL 

body  at  ambient  temperature. 

Equation  1-1  is  then  written  in  terms  of  a  "mult"  and  an  "add" 
factor  as 


L*  -  VT  •  MULT  +  ADD 


(1-4) 


where 


MULT  -  W 

ADD  -  L®B-  WVa 


W  -  (l£  -  lad)/(v£  -  vA) 


(1-5) 
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The  computer  program  used  for  calibration  initially  averages  the 
integer  values  (bin  values)  of  the  calibration  sources  for  all  scan  lines 
in  the  image  and  outputs  a  mean  and  a  standard  deviation  for  them.  If 
the  standard  deviations  are  small  compared  to  the  mean  values  the  program 
allows  for  an  "average"  calibration  in  which  the  mean  bin  values  of 
and  VA  are  used  in  Equations  1-4  and  1-5.  In  this  case  the  integer 
data  values  on  the  calibrated  tape  are  the  same  as  those  in  the  original 
tape;  only  "mult"  and  "add"  factors  are  evaluated  by  the  calibration 
procedure.  However,  if  the  standard  deviations  are  not  small  compared 
to  the  means,  indicating  a  drift  in  the  system,  a  "line-by-line"  cali¬ 
bration  is  required.  For  this  type  of  calibration  values  of  and 
are  taken  from  the  calibration  sources  at  the  end  of  each  scan  line  and 
these  values  used  to  determine  "mult"  and  "add"  factors  from  Equation 
1-5.  These  factors  are  then  used  to  modify  the  integer  data  values  of 
a  given  scan  line  so  that  only  one  "mult"  and  "add"  factor  is  required 
for  the  total  image.  The  single  "mult"  and  "add"  factors  used  are  those 
determined  from  the  mean  V^,  values  and  the  modified  integer  data 
value  is  taken  as 


VT  •  MULT  +  ADD  -  ADD 

WET 


(1-6) 


where  the  average  factors  are  MULT  and  ADD,  those  determined  for  each 

I 

line  are  MULT  and  ADD,  and  VT  is  the  modified  integer  data  value. 

The  long  wavelength  IR  channels  (X  >  3.0  pm)  are  calibrated  in 
temperature  and  require  a  somewhat  more  involved  calibration  procedure. 
Since  the  detector  output  is  proportional  to  radiance  and  not  temperature, 
the  first  step  in  the  process  is  to  relate  bin  value  (integer  data  value) 
to  radiance,  as  in  the  reflective  IR  channels,  using  the  "cold"  plate 
as  a  dark  level  reference,  the  "hot"  plate  as  a  radiance  standard,  and 
the  "ambient"  plate  to  correct  for  ambient  radiation  levels.  If  a 
linear  relationship  is  again  assumed  between  detector  output  and  radiance 
the  apparent  target  radiance  is  given  by  an  equation  identical  to  1-1 
except  that  the  ambient  radiance  (L^)  and  bin  value  (V  )  are  replaced 


( 
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by  those  for  the  "cold"  plate  (L^,  V^)  and  the  lamp  radiance  (L^)  and  bin 

value  (V»)  are  replaced  by  those  for  the  "hot"  plate  (L  ,  V  ) : 

H  H 


/vi  ■  vc\ 

Lc  +  (VH  -  VC  )  <Lh  "  Lc> 


(1-7) 


In  this  case  however  the  "hot"  plate  radiance  is  given  by 


T  T  /I  \  T 

Si  “  EHLH  +  (1  “  LA 


(1-8) 


where  e  is  the  emissivity  of  the  "hot"  plate,  L„  the  radiance  of  a 

"  **  BB 

blackbody  at  the  same  temperature  as  the  "hot"  plate,  and  the  am¬ 
bient  blackbody  radiance.  Similarly,  the  "cold"  plate  radiance  is  given 
by 


t  ,  /i  .  T BB 

ecLc  +  (1  "  ec)  la 


(1-9) 


In  both  of  these  expressions  the  first  term  is  the  radiance  emitted  by 
the  plate  while  the  second  term  is  the  ambient  radiance  reflected  by  the 
plate. 

Substitution  of  Equations  1-8  and  1-9  into  1-7  and  taking  the 
plate  emissivities  equal  gives 

,  a  f  BB  .  .BBl  ,  /VS  ~  Vc\  f  ,,  BB  _  BB/l  ,T-m\ 

LT  "  [£LC  (1  ”  J  LA  J  ^Vh  -  Vc  J  [e(LH  “  LC  (I”10) 

This  equation  could  be  used  for  calibration  but  the  ambient  radiance 
BB 

would  have  to  be  calculated  from  ambient  temperature  measurements. 
This  may  be  avoided  by  again  assuming  linearity  of  the  detectors  and 
taking 


,BB 

la  *  LC  + 


ft3) 


<lh  -  Lc> 


(1-11) 


as  for  L*  in  Equation  1-7.  Substituting  Equations  1-8  and  1-9  into  this 
expression  it  is  found  that  many  of  the  terms  cancel  and  the  ambient 
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radiance  is  given  by 


LBB  -  LBB  + 

la  lc 


V  -  V  \ 

_A _ C\  .  BB  BB. 

c> 


(1-12) 


where  all  radiances  are  now  blackbody  radiances  evaluated  at  the  plate 
temperatures.  Making  a  final  substitution  of  Equation  1-12  into 
Equation  1-10  and  rearranging  terms  gives 


Lt  -  L*B  +  W  [(1  -  e)VA  -  Vc]  +  eWV1 


(1-13) 


where 


(1-14) 


Equation  1-13  then  relates  the  data  bin  values  (V^,)  to  radiance  and  a 
second  transformation  is  required  to  convert  these  values  to  temperature. 

The  computations  required  to  convert  radiance  to  temperature  have 
been  minimized  by  noting  that  the  radiance  values  are  bounded;  the 
minimum  value  being  given  by  =  0  and  the  maximum  by  =  255.  Using 
this  fact  a  table  of  temperature  versus  radiance  is  constructed,  as 
shown  in  Table  1-1,  by  integvating  the  Planck  function  over  the  bandpass 
of  the  channel  to  be  calibrated  choosing  a  range  of  temperatures  suf¬ 
ficient  to  cover  the  entire  radiance  range.  The  temperature  increments 
in  this  table  are  taken  small  enough  (~1°C)  so  that  interpolation  may 
be  used  to  form  an  array  of  temperature  versus  data  bin  value  T(V^)  for 
all  possible  bin  values.  It  is  this  array  which  is  then  used  to  directly 
convert  the  uncalibrated  data  values  to  temperature. 

Unlike  the  radiance  data,  the  integer  values  appearing  in  the 
calibrated  images  for  the  thermal  channels  are  never  the  same  as  those 
of  the  original  image  since  the  calibrated  values  are  now  proportional 
to  temperature  not  radiance.  This  new  set  of  integer  data  values  is 
scaled  from  0  to  255  using 
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-la.  TABLE  OF  BLACXBODY  RADIANCE  vs  TEMPERATURE  IN  THE  4.5  -  5.5  urn  BAND 
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TABLE  I-lb.  TABLE  OF  BLACKBODY  RADIANCE  VS  TEMPERATURE  IN  THE  9.3  -  11.7  pm 
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TABLE  I-lc.  TABLE  OF  BLACKBODY  RADIANCE  vs  TEMPERATURE  IN  THE  8.0-13.5  mid  BAND 


T (V  )  -  1(0) 

N  =  1(255)  -  T (0)  255  (I~ 

where  1  is  the  temperature  versus  bin  value  array  with  1(VT)  equal  to 
the  temperature  corresponding  to  bin  value  V^.  The  "mult"  and  "add" 
factors  for  the  thermal  channels  are  then  taken  as 


MUL1 


1(255)  -  1(0) 
255 


ADD  -  1(0) 


(1-16) 


Ihe  calibration  program  allows  for  "average"  or  "line-by-line" 
calibration  of  the  thermal  channels  as  it  did  for  the  reflective  IR 
channels.  Because  of  the  amount  of  computation  required  to  set  up  the 
temperature  versus  bin  value  array,  only  one  such  array  is  generated 
using  the  scene  average  values  of  V  ,  V  ,  and  V  in  Equation  1-13.  If 

LA  rl 

an  "average"  calibration  is  used  this  array  directly  relates  data  bin 
values  (V^)  to  temperature  and  Equation  1-15  is  used  directly.  However, 
if  a  "line-by-line"  calibration  is  required  the  data  bin  values  must  be 
modified  so  that  1-13,  using  scene  average  values  of  V  ,  V  ,  and  V  , 
gives  the  correct  radiance.  Ihe  modification  of  the  data  bin  values 
is  done  in  the  same  way  as  for  the  reflective  IR  channels  using  Equation 
1-6.  In  this  case  the  "mult"  factors  are.  given  by  eW  or  eW  and  the 
"add"  factors  by  the  first  two  terms  of  Equation  1-13  or 

ADD  =  LBB  +  W  [(1  -  e)  VA  -  Vc] 

ADD  -  LBB  +  W  [(1  -  z)  -  V^] 

lo  avoid  round-off  errors  the  values  determined  from  1-6  are  not  rounded 
to  the  nearest  integer  but  taken  as  floating  point  numbers  and  inter¬ 
polation  in  the  I(VT)  array  is  used  to  determine  the  temperature.  Ihis 
temperature  is  then  used  in  place  of  1(VT)  in  Equation  1-15  to  evaluate 
the  calibrated  integer  value. 


In  a  typical  data  set  both  thermal  and  reflective  IR  channels  are 
present  so  that  both  of  the  above  procedures  are  used  in  generating  a 
calibrated  image  tape. 


CORRECTIONS  FOR  OVERLAP  AND  UNEQUAL  FIELDS-OF-VIEW 

Since  the  detectors  used  in  the  various  channels  (wavelength 

regions)  vary  in  size  and  consequently  in  instantaneous  field -of -view, 
direct  comparisons  of  the  channels  is  not  possible  unless  the  fields- 
of-view  are  normalized  in  some  fashion.  A  second,  although  less  critical, 
problem  is  that  depending  on  the  aircraft  altitude,  scan  lines  and/or 
pixels  may  be  highly  overlapped  resulting  in  oversampling  of  the  scene. 
Both  of  these  problems  were  corrected  at  the  time  the  data  was  calibrated 
by  averaging  the  scan  lines  (and/or  pixels)  of  the  smaller  f ield-of-view 
channels  and  generating  an  equivalent  scan  line  (and/or  pixel)  with  a 
f ield-of-view  equal  to  the  largest  f ield-of-view.  This  procedure  must 
be  carried  out  separately  for  scan  lines  and  pixels  within  a  scan  line 
since  the  fields-of-view  and  the  data  taking  rates  are  different  in 
these  two  dimensions.  In  the  following  discussion  only  scan  line, 
averaging  will  be  discussed  since  pixel  averaging  is  completely 
analogous . 

As  outlined  in  Section  3.1,  the  averaging  procedure  consists  of 
setting  up  a' set  of  contiguous  output  lines  with  a  f ield-of-view  equal 
to  the  largest  f ield-of-view  (D)  and  then  combining  the  smaller  field- 
of-view  (d)  scan  lines  to  generate  equivalent  lines  with  fields-of-view  D. 
As  shown  in  Figure  1-1,  the  combining  of  scan  lines  is  accomplished 
using  a  weighted  average  given  by 

Lf  v 

V  _  (I_17) 


where  V^  is  the  k  output  line  data  value  for  a  given  pixel,  the 
input  line  data  value  for  the  same  pixel,  and  F^  ^  the  overlap 
between  input  line  t  and  output  line  k.  Since  the  overlap  will  be  non¬ 
zero  only  for  a  small  number  of  lines  about  the  output  line,  F^  ^  need 
only  be  calculated  for  a  few  values  of  l .  A  general  expression  for  the 
overlap  was  determined  by  considering  the  displacements  of  the  scan  lines 
from  the  beginning  of  the  scene.  The  displacement  to  the  bottom  of 
input  scan  line  l  is  given  by 
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INPUT  OUTPUT 


FIGURE  1-1  SCHEMATIC  REPRESENTATION  OF  THE  PROCEDURES  USED  FOR  SCAN  LINE  AVERAGING 


(1-18) 


U  -  1)  VaT  + 


(D  +  d) 
2 


where  is  the  aircraft  ground  velocity  and  t  the  scan  period.  The 
second  term  in  this  equation  is  included  since  the  smaller  field-of-view 
(d)  is  concentric  with  D  in  scan  line  1  and  consequently  has  an  initial 
displacement  (to  the  bottom  of  the  scan  line)  given  by 


(1-19) 


As  shown  in  Figure  1-2  the  overlap  between  input  line  t  and  output  line 
k  is  then  given  by 


Y 


a  -  Dvax  + 


-  (k  -  1)D 


(1-20) 


The  absolute  value  appears  in  this  equation  to  account  for  cases  where  the 
bottom  of  the  input  line  has  a  greater  displacement  than  that  of  the 
outj-  ».  ine.  It  may  be  seen  that  the  relationship  between  k  and  l  is 

V  T 

(k  -  1)  =  Integer  value  of  [(£  -  1)  ]  (1-22) 


These  overlap  factors  are  then  used  in  Equation  1-17  to  generate  the 
equivalent  output  lines.  Since  in  general  the  scan  lines  for  the  largest 
field-of-view  channel  are  themselves  overlapped,  this  procedure  is  used 
to  average  these  lines  as  well  so  that  a  contiguous  set  of  output  lines 
is  obtained  in  all  channels.  In  this  case,  the  same  equations  are  used 
but  d  is  equal  to  D. 
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FIGURE  1-2  SCHEMATIC  OF  DISPLACEMENTS  USED  TO  DERIVE  THE  OVERLAP  FACTORS  F 
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APPENDIX  II 

COMPLETE  COMPILATION  OF  STATISTICAL  DATA 

DISCUSSION  OF  FORMATS,  SPECIFIC  PROCEDURES,  ANOMALIES  AND  LIMITATIONS 
This  appendix  presents  the  statistics  generated  for  all  of  the 
chosen  scenes.  It  is  broken  down  into  two  major  sections,  the  first 
of  which  presents  the  mean  values,  standard  deviations,  spectral 
correlations,  and  histograms  for  each  of  the  scenes,  while  the  second 
section  presents  the  area/intensity  statistics  and  Wiener  spectra 
computed  for  Flint-1  and  Mill  Creek.  It  is  important  to  emphasize  that 
the  data  presented  here  are  only  descriptive  of  the  specific  scenes 
for  the  specific  conditions  under  wl^.ch  the  imagery  was  collected  and 
that  generalizations  to  other  background  scenes,  other  background  types 
or  other  conditions  cannot  be  made  without  further  data  analysis  and 
modeling.  The  primary  purpose  of  this  presentation  is  to  familiarize 
the  potential  user  of  the  data  with  what  is  and  what  can  be  made 
available  to  assist  in  systems  analysis. 

In  the  first  section  of  this  appendix,  the  statistics  for  each 
scene  are  preceded  by  (1)  maps  showing  the  regions  covered,  (2)  line- 
scan  images  for  each  of  the  wavelength  bands  generated  from  the  scanner 
data,  and  (3)  a  visual  image  showing  the  subareas  defined  for  statistics 
generation.  The  histograms  for  the  subareas,  total  area,  and  calibra¬ 
tion  sources  of  each  scene  are  presented  in  graphical  form  but  tabu¬ 
lations  of  the  data  used  to  generate  these  plots  are  available  at 
ERIM  for  quantitative  analysis.  Histograms  of  the  calibration  sources 
have  been  included  for  two  reasons: 

(1)  To  indicate  the  range  of  possible  data  values  so 
saturation  problems  may  be  identified,  and 

(2)  To  indicate  the  distribution  of  calibration  source 
values  and  hence  the  possible  uncertainty  associated 
with  the  radiance  or  temperature  calibrations. 

Histograms  of  the  calibration  sources  are  not  presented  for  the  9.3- 
11.7  ym  Baltimore  data,  the  1.0-1. 4  and  2. 0-2. 6  urn  Black  Hills-1  data, 
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and  the  1.0-1. 4  and  2. 0-2. 6  pm  Mono  Lake  data.  Only  the  position  of 
the  calibration  signal  used  is  shown.  The  Baltimore  data  were  processed 
before  the  calibration  histogram  capability  was  developed,  and  the 
Black  Hills-1  and  Mono  Lake  calibration  histograms  are  not  shown  because 
the  calibration  was  inferred  from  other  data  as  discussed  later. 

Each  of  the  calibration  source  histograms  is  presented  with  dual 
abscissa  units:  the  recorded  integer  data  value  (bin  number)  and  the 
calibrated  radiance  or  temperature.  Since  the  data  is  collected  as 
positive  8-bit  integers,  this  data  is  constrained  to  the  range  0  to  255 
with  negative  data  points  being  given  the  value  zero  and  points  with 
values  greater  than  255  being  given  the  value  255.  Consequently,  if 
the  wings  of  a  given  data  distribution  show  the  presence  of  scene 
elements  with  calibrated  data  valyes  corresponding  to  a  bin  value  of 
zero  or  255,  it  must  be  assumed  that  these  points  are  saturated  and 
that  the  actual  scene  contained  values  outside  the  0  -  255  range. 

Cases  where  255  saturation  are  in  evidence  are  in  the  Flint-2  2. 0-2. 6  pm 
channel,  the  Mill  Creek  1.5-1. 8  and  2. 0-2. 6  pm  channels,  the  Black  Hills-1 
2. 0-2. 6  pm  channel,  and  the  Black  Hills-2  1. 5-1.8  pm  channel.  Zero 
saturations  were  found  in  the  Baltimore  9.3-11.7  pm  channel. 

The  histograms  of  the  calibration  sources  cannot  be  used  to  directly 
infer  the  influence  of  system  noise  on  the  distribution  of  scene  data 
values.  This  is  the  case  since  in  many  instances  the  spread  of  cali¬ 
bration  source  values  is  due  to  a  cyclical  variation  or  a  drift  of  the 
source  throughout  the  scene  which  is  corrected  for  using  the  line-by¬ 
line  calibration  procedures  discussed  in  Section  3.  No  attempt  has 
been  made  at  the  present  time  to  determine  the  degree  to  which  these 
line-by-line  procedures  completely  removed  these  cyclical  or  drift 
variations  during  calibration. 

A  few  comments  regarding  the  data  appearing  in  the  first  section 
of  this  appendix  are  in  order: 

(1)  Each  of  the  histograms  for  the  subareas  and  total  area  are 
independently  normalized  by  assigning  a  relative  frequency 
of  unity  to  the  data  value  appearing  most  frequently  in  the 
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scene.  As  a  result,  if  255  saturations  occur  and  a  large 
enough  number  of  points  have  this  value,  the  histogram  may 
never  reach  unity  except  for  the  spike  appearing  at  255. 

This  is  in  evidence  in  the  histogram  of  the  2. 0-2. 6  vim  channel 
of  Flint-2  data. 

(2)  In  addition  to  the  digital  saturations  discussed  earlier, 
some  of  the  data  histograms  show  an  anomalously  large  number 
of  data  points  at  the  high  end  of  the  distribution  which  are 
not  digitally  saturated.  In  the  cases  where  this  has  been 
identified,  it  was  found  to  be  due  to  analog  saturation  in 
the  aircraft  recording  system,  and  hence  is  not  indicative 
of  the  true  scene  distribution.  Cases  where  this  has  been 
observed  are  the  Flint-1  9.3-11.7  um  channel  around  data 
value  308K,  the  Black  Hills-1  4. 5-5. 5  urn  channel  near  data 
value  307K,  and  the  Black  Hills-2  2. 0-2. 6  pm  channel  near 
data  value  280  yw/ cm2 • sr • ym) . 

(3)  In  the  9.3-11.7  ym  channels  of  Flint-1,  Flint-2,  and  Mill 
Creek,  as  well  as  the  11.3-13.5,  8.0-10.9,  and  9.4-12.1  ym 
channels  of  Pisgah  Crater,  some  spikes  appear  in  the  data 
histograms.  These  are  the  results  of  the  analog-to-digital 
conversion  process  in  which  certain  data  values  are  digitized 
as  the  next  higher  or  lower  integer  value.  This  digiti¬ 
zation  error  is  not  important  as  it  is  an  error  in  magnitude 
of  only  a  few  tenths  degree  Kelvin.  These  are  not  seen  in 
the  imagery  since  these  data  values  appear  randomly  through¬ 
out  the  scene  but  they  do  appear  in  the  histograms. 

"ft 

As  discussed  elsewhere  in  this  report  ,  the  computer  code  developed 
for  data  calibration  was  equipped  to  do  three  different  types  of 
calibration: 


I 


*See  Section  3  and  Appendix  I. 
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(1)  An  average  calibration  in  which  the  average  bin  values 
of  the  calibration  sources  for  the  whole  scene  are  used; 

(2)  A  line-by-line  dark  level  correction  in  which  the  dark 
level  value  of  each  scan  line  is  used  in  conjunction 

with  average  values  for  the  other  calibration  sources;  and 

(3)  A  line-by-line  calibration  in  which  the  bin  values  of 
all  calibration  sources  taken  from  the  individual  scan 
lines  are  used. 

In  general,  an  average  calibration  was  used  when  the  calibration 
sources  were  reasonably  constant  throughout  the  scene  so  that  their 
average  bin  values  had  small  standard  deviations,  while  a  line-by-line 
procedure  was  used  when  the  calibration  sources  showed  cyclical  varia¬ 
tions  or  drifts.  Of  the  two  possible  line-by-line  procedures,  a  dark 
level  correction  (2)  was  used  if  only  the  dark  level  varied,  while  a 
complete  line-by-line  calibration  (3)  was  used  if  any  or  all  of  the 
other  calibration  sources  showed  variation.  However,  in  its  present 
form,  the  computer  code  is  constrained  to  use  the  same  type  of  cali¬ 
bration  on  all  channels  calibrated  in  radiance  and  a  separate  type  for 
all  channels  calibrated  in  temperature.  Consequently,  if  any  of  the 
radiance  (or  temperature)  channels  required  a  line-by-line  calibration, 
the  remainder  of  the  radiance  (or  temperature)  channels  had  to  be 
similarly  calibrated. 

To  assist  in  determining  the  calibration  procedure  required  for 
each  scene,  three  types  of  data  displays  have  been  used: 

(1)  Histograms  of  the  calibration  sources  for  the  total  scene 

to  determine  mean  values  and  standard  deviations  (Figure  I I— 1 ) ; 

(2)  "A"-scope  traces  of  several  full  scan  lines  to  determine 
line-to-line  variations  and  to  see  if  variations  in  the  video 
portion  of  the  scan  followed  those  in  the  calibration  sources 
(Figure  II-2);  and, 

(3)  Time  history  traces  of  the  calibration  sources  to  observe 
source  drifts  or  cyclical  variations  (Figure  II-3). 
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FIGURE  I 1-3.  EXAMPLES  OF  TIME  HISTORY  TRACES  OF  THE  COLD,  AMBIENT,  AND  HOT  PLATE 

DATA  VALUES  FOR  A  THERMAL  IR  CHANNEL 


Table  II-l  outlines  the  calibration  procedures  used  for  each  of  the 
data  sets  along  with  comments  explaining  the  choice  of  procedures  used. 
The  Black  Hills-1  and  Mono  Lake  data  were  the  only  scenes  requiring 
special  calibration  procedures. 

The  data  for  Black  Hills-1  was  gathered  using  the  M5  scanner  with 
the  thermal  calibration  plates  in  part  of  the  f ield-of-view  which,  at 
the  time  this  data  was  collected,  did  not  also  have  radiance  calibration 
sources.  Hence,  the  Black  Hills-1  radiance  channels  could  not  be 
calibrated  directly.  However,  the  data  for  Black  Hills-2  was  gathered 
at  the  same  time  as  Black  Hills-1  using  the  other  M5  scanner  which  did 
have  the  required  calibration  sources.  Since  both  of  these  data  sets 
had  1.0-1. 4  pm  and  2. 0-2. 6  pm  channels,  the  Black  Hills-1  data  in  these 
channels  were  calibrated  by  equating  the  average  radiance  in  a  given 
region  of  the  scene  with  that  of  Black  Hills-2  for  the  same  region  of 
the  scene . 

The  Mono  Lake  data  did  contain  the  required  radiance  calibration 
sources  but  the  time  history  traces  showed  that  these  sources  were  not 
stable  but  varied  randomly  throughout  the  scene.  Since  the  "A"-scope 
traces  indicated  that  the  video  was  stable  and  did  not  follow  the 
calibration  source  fluctuations,  it  was  decided  that  the  sources  could 
not  be  used  for  reliable  calibration.  Fortunately,  the  Mono  Lake  and 
Black  Hills  data  sets  were  gathered  with  the  same  scanner  systems  but 
with  amplifier  gains  differing  by  a  factor  of  two.  Assuming  that  the 
system  response  is  reasonably  constant,  the  1.0-1. 4  and  2. 0-2. 6  pm 
Mono  Lake  data  were  calibrated  using  the  response  values  obtained  from 
the  Black  Hills  data  sets  corrected  by  the  difference  in  amplifier 
gains.  Because  the  assumption  of  system  stability  in  the  radiance 
channels  could  not  be  tested  directly,  the  possible  errors  in  this 
calibration  could  not  be  estimated.  The  thermal  channels  (4. 5-5. 6  pm 
and  8.0-13.5  pm)  of  Mono  Lake  did  not  have  this  problem  and  were 
calibrated  normally. 

The  second  section  of  this  appendix  presents  one-dimensional 
Wiener  spectra  for  the  Flint-1  and  Mill  Creek  scenes,  a  representative 
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TABLE  I I— 1 .  CALIBRATION  PROCEDURES  USED  FOR  THE  VARIOUS  BACKGROUND  SCENES 
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two-dimensional  Wiener  spectrum  for  the  9.3-11.7  pm  band  of  Flint-1, 
and  area/ intensity  statistics  for  the  1.0-1. 4  pm  and  9.3-11.7  pm  bands 
of  Flint-1  and  Mill  Creek. 

Included  with  the  one-dimensional  Wiener  spectra  for  Flint-1  are 
spectra  of  the  dark  level  (ambient  or  cold  calibration  source)  since 
these  plots  indicate  that  the  feature  at  0.5  cycles/meter  is  not  a  scene 
feature  but  a  spatial  frequency  generated  by  the  scanner  itself.  The 
area/ intensity  statistics  are  presented  in  two  forms:  plots  of  the 
equivalent  elliptical  areas  for  each  threshold  and  tabulations  of 
these  areas  as  sorted  by  geometric  area,  perimeter,  and  shape  factor. 
These  latter  tabulations  are  broken  down  into  coarse  increments,  only 
to  cover  a  large  enough  range  to  be  generally  applicable;  these  incre¬ 
ments  are  not  indicative  of  the  resolution  of  the  system  and  tabulations 
with  much  finer  increments  may  be  generated  as  required. 

Calibrated  data  tapes  for  each  of  the  scenes  presented  have  been 
prepared  and  are  available  for  additional  processing  as  the  need  for 
additional  statistics  becomes  apparent. 
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MEAN  VALUES,  STANDARD  DEVIATIONS,  SPECTRAL 
CORRELATIONS,  AND  HISTOGRAMS  FOR  THE  SUB- 
AREAS  AND  TOTAL  AREAS  OF  THE  CHOSEN  BACK¬ 
GROUND  SCENES 


I  mCEDUO  BIOS  BUMrnNOT  71LMED 


4r 


* 

FLINT-1 


Scene  Type  Residential 

Date  of  Flight  18  September  1971 

Time  of  Flight  1132  -  1133 

Altitude  (Ft)  1000 

No.  of  Sub-Areas  6 

No.  of  Data  Points  514,065 


Channels 

2 

3 

4 

5 

Wavelength  (um) 

1.0-1. 4 

1.5-1. 8 

2. 0-2. 6 

9.3-11. 

Resolution  (mr) 

In-Track 

5.0 

5.0 

5.0 

5.0 

Cross-Track 

2.5 

2.5 

2.5 

2.9 

Nadir  Pixel  Dimensions 
In-Track 

(m) 

1.524 

1.524 

1.524 

1.524 

Cross-Track 

0.762 

0.762 

0.762 

0.884 

Nadir  Ground  Sample 
Distance  (m) 

In-Track 

1.524 

1.524 

1.524 

1.524 

Cross-Track 

0.762 

0.762 

0.762 

0.762 

Line  Averaging  used  on  ALL  channels. 


*These  data  were  obtained  with  the  M-7  scanner.  All  data  are  in 
spatial  registration. 
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Pixel  1 


Pixel  161 


Line  403 


SUB-AREAS  DEFINED  FOR  STATISTICS  GENERATION  IN  THE  FLINT-1  IMAGE. 
Uneven  areas  were  chosen  so  that  Areas  2  and  5  covered  the  +20° 
range  suitable  for  correlation.  Approximate  scene  dimensions  are 
3985  ft  (1215  m)  by  1612  ft  (491  m) .  Each  sub-area  as  well  as  the 
total  area  have  been  histogrammed.  Histogram  plots  and  their 
respective  sub-areas  are  identified  with  the  following  key: 


FLINT-1 
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IUTAL  PIS.  M82P.  63828.  63B2B.  63828 


*u-*m  38*  ENVIRONMENTAL  RESEARCH  INST  of  MICHIGAN  ANN  ARBOR  IN— ETC  F/6  17/5 
STATISTICAL  ANALYSIS  OF  TERRAIN  BACKGROUND  MEASUREMENTS  DATA* (U) 

MAR  77  R  SPELLICT  »  J  BEARD  »  J  R  MAXWELL  N00123-76-C-0708 

UNCLASSIFIED  ERIM-120500-12-F 
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fWEflt  FLINT  1 
LAMBDA*  9.3  TO  11.7  mM 


CflLJB.PLflTES 


FLINT-2* 

Industrial 
18  Sep  1971 
1155  -  1157 
1000 
6 

572,115 


Channels 

2 

3 

4 

5 

Wavelength  (pm) 

1.0-1. 4 

1.5-1. 8 

2 

.0-2.6 

9.3-11.7 

Resolution  (mr) 

In-Track 

5.0 

5.0 

5.0 

5.0 

Cross-Track 

2.5 

2.5 

2.5 

2.9 

Nadir  Pixel  Dimensions  (m) 

In-Track 

1.524 

1.524 

1.524 

1.524 

Cross-Track 

0.762 

0.762 

0.762 

0.884 

Nadir  Ground  Sample 

Distance  (m) 

In-Track 

1.524 

1.524 

1.524 

1.524 

Cross-Track 

0.762 

0.762 

0.762 

0.762 

Line  averaging  used 

on  all 

channels 

of 

data 

Scene  Type 
Date  of  Flight 
Time  of  Flight 
Altitude  (Ft.) 

No.  of  Sub-Areas 
No.  of  Data  Points 


*  These  data  were  obtained  with  the  M-7  scanner.  All  data  are 
on  spatial  registration. 
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LINE  SCAN  DIAGES  PRODUCED  FROM  THE  VARIOUS  INFRARED  CHANNELS  OF  FLINT-2 
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SUB-AREAS  DEFINED  FOR  STATISTICS  GENERATION  IN  THE  FLINT-2  IMAGE. 
Uneven  areas  were  chosen  so  that  Areas  2  and  5  covered  the  ±20° 
range  suitable  for  correlation.  Approximate  scene  dimensions  are 
4435  ft  (1352  m)  by  1612  ft  (491  m).  Each  sub-area  as  well  as  the 
total  area  have  been  hlstogrammed.  Histogram  plots  and  their 
respective  sub-areas  are  identified  with  the  following  key: 

O  Sub-area  1  +  Sub-area  4 


O  Sub-area  2 


x  Sub-area  5 


A  Sub-area  3 


O Sub-area  6 


11-48 


+  + 

Ll  ll 

r-  x 

a.  r* 

a  ** 

•  • 

rn  — 


c\l 

H 

2 

HH 

►J 


W 

o 

a 


g 

H 


f\. 


G 

f'j 

— 

c 

4- 

*• 

+ 

• 

OJ 

u. 

«— 

G 

K*. 

a 

r- 

a 

ct 

G 

K* 

r^ 

G 

fO 

• 

• 

O 

o 

0.- 

r* 

• 

• 

*— 

G 

1 

c 

G 

g 

K* 

X 

rvc 

(\j 

c 

>C 

-* 

G 

G 

• 

• 

• 

4- 

+ 

— 

c 

c 

li. 

Ll 

1 

h* 

X 

r\. 

a 

G 

c 

O' 

a 

lT 

o 

X 

Kl 

X 

G 

G 

G 

G 

• 

• 

• 

• 

• 

• 

rr 

G 

o 

O 

IIITAL  PTS.  S7?U5  STPIfS  S7P  I  1 S  S  /  P  1  I  r> 

HIM!  2 


11-50 


0.00  500.00  1000.00  1500.00  2000.00  2500.00  3000.00 

wW/CM2  STR  juM 


0.00  500.00  1000.00  1500.00  2000.00  2500 

mW/CM2  STR  wM 


CflLIB. PLATES 


nW/CM2  STR  mM 


11-54 


11-56 


11-57 


AREA*  FLINT  2 


flREflj  FLINT  2 


PRECEDING  PAQ 1  BLANK. NOT  FTl>,gD 


T 

BALTIMORE* 


Scene  Type 
Date  of  Flight 
Time  of  Flight 
Altitude  (Ft) 

No.  of  Sub-Areas 
No.  of  Data  Points 


Residential 
11  May  1972 
1137  -  1139 
2500 
6 

269,610  for  channels  2,4 
322,500  for  channel  5 


Channels 

2 

4 

5 

Wavelength  (pm) 

1.0-1. 4 

2. 0-2. 6 

9.3-11. 

Resolution  (mr) 

In-Track 

5.0 

5.0 

5.0 

Cross-Track 

2.5 

2.5 

2.9 

Nadir  Pixel  Dimensions 
In-Track 

(m) 

3.810 

3.810 

3.810 

Cross-Track 

1.905 

1.905 

2.210 

Nadir  Ground  Sample 

Distance  (m) 

In- Track 

3.810 

3.810 

3.810 

Cross-Track 

1.905 

1.905 

1.905 

Line  Averaging  used  for  ALL  channels. 


*These  data  were  obtained  with  the  M-7  scanner.  The  1.0-1. 4  and 
2.0-2. 6  pm  data  are  in  spatial  registration,  but  the  9.3-11.7  pm 
data  were  processed  separately  and  are  not  in  spatial  registration 
with  the  1.0-1. 4  and  2. 0-2. 6  pm  data.  Hence,  spectral  correlation 
coefficients  have  not  been  determined  between  the  9.3-11.7  pm  data 
and  either  the  1.0-1. 4  pm  or  2. 0-2. 6  pm  data. 
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LINE  SCAN  IIIAGES  PRODUCED  FROM  THE  VARIOUS 
INFRARED  CHANNELS  OF  BALTIMORE 
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Line  427  -  Channels  2,  4 
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SUB-AREAS  DEFINED  FOR  STATISTICS  GENERATION  IN  THE  BALTIMORE  IMAGE. 
Uneven  areas  were  chosen  so  that  Areas  2  and  5  covered  the  ±20° 
range  suitable  for  correlation.  Approximate  scene  dimensions  are 
6250  ft  (1905  m)  by  4031  ft  (1229  m)  for  channel  5  and  5225  ft 
(1592  m)  by  4031  ft  (1229  m)  for  channels  2  and  4.  Each  sub-area 
as  well  as  the  total  area  have  been  histogrammed.  Histogram  plots 
and  their  respective  sub-areas  are  identified  with  the  following  key 
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MILL  CREEK 


Scene  Type 
Date  of  Flight 
Time  of  Flight 
Altitude  (Ft) 

No.  of  Sub-Areas 
No.  of  Data  Points 

Channels  2 


Wavelength  (  m)  1.0-1. 4 

Resolution  (mr) 

In-Track  5.0 

Cross-Track  2.5 

Nadir  Pixel  Dimensions  (m) 

In-Track  4.572 

Cross  Track  2.286 

Nadir  Ground  Sample 
Distance  (m) 

In-Track  4.572 

Cross-Track  2.286 


Mountainous 
30  June  1972 
0733  -  0736 
3000 
6 

171,570 


3 

4 

5 

5-1.8 

2. 0-2. 6 

9.3-11. 

5.0 

5.0 

5.0 

2.5 

2.5 

2.9 

4.572 

4.572 

4.572 

2.286 

2.286 

2.652 

4.572 

4.572 

4.572 

2.286 

2.286 

2.286 

Line  Averaging  used  for  ALL  channels. 


*These  data  were  obtained  with  the  M-7  scanner.  All  data  are 
in  spatial  registration 
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MAP  OF  WICHITA  MOUNTAINS  AREA  SHOWING  THE  REGION  COVERED  IN  THE  MILL  CREEK  RUN  ABOUT  FLIGHT  LINE  6 


SUB-AREAS  DEFINED  FOR  STATISTICS  GENERATION  IN  THE  MILL  CREEK  IMAGE. 
Uneven  areas  were  chosen  so  that  Areas  2  and  5  covered  the  ±20° 
range  suitable  for  correlation.  Approximate  scene  dimensions  are 
3990  ft  (1216  m)  in-track  by  4837  ft  (1474  m)  cross-track.  Each 
sub-area  as  well  as  the  total  area  have  been  histogrammed.  Histogram 
plots  and  their  respective  sub-areas  are  identified  with  the  following 
key: 
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AREAt  HILL  CHEEK 

LAM0OA-  1.0  TO  1.4  pH  MEAN  -  112.2! 

TOTAL  AREA  ST. DEV.-  62.2 
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BLACK  HILLS-1* 


Scene  Type 
Data  of  Flight 
Time  of  Flight 
Altitude  (Ft) 

No.  of  Sub-Areas 
No.  of  Data  Points 


Forested  Mountains 
22  July  1969 
1340  -  1342 
1500 
4 

224,130  for  channels  2,  4,  and  5 
425,630  for  channel  12 


Channels 

2 

4 

5 

12 

Wavelength  (um) 

1.0-1. 4 

2.0-2. 6 

4. 5-5. 5 

8.0-13. 

Resolution  (mr) 

In-Track 

6.6 

6.6 

6.6 

3.5 

Cross-Track 

3.5 

3.5 

3.5 

3.5 

Nadir  Pixel  Dimension 
In-Track 

(m) 

3.017 

3.017 

3.017 

1.600 

Cross-Track 

1.600 

1.600 

1.600 

1.600 

Nadir  Ground  Sample 

Distance  (m) 

In-Track 

3.017 

3.017 

3.017 

1.600 

Cross-Track 

1.143 

1.143 

1.143 

1.143 

Line  Averaging  used  for  channels  2,  4,  and  5  and  for  channel  12 
independently. 


*The  Black  Hills-1  data  were  collected  with  an  M-5  scanner  with 
thermal  calibration  plates  in  part  of  the  scanner  f ield-of-view. 

The  8.0-13.5  um  detector  and  the  1.0-1. 4,  2. 0-2. 6,  4.5-5.5um 
detectors  were  on  opposite  ends  of  the  scanner  and  are  not  in  spatial 
registration.  Hence,  spectral  correlation  coefficients  have  not  been 
determined  between  the  8.5-13.5  um  data  and  the  1.0-1. 4,  2. 0-2. 6,  or 
4.5-5. 5  um  data. 
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MAP  OF  BLACK  HILLS  1 


LINE  SCAN  IMAGES  PRODUCED  FROM  THE  VARIOUS 
INFRARED  CHANNELS  OF  BLACK  HILLS-1 


Line  1  Line  344 
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Line  1034  Line  1374  -  channel  12 

Line  552  Line  732  -  channels  2,4,5 


SUB-AREAS  DEFINED  FOR  STATISTICS  GENERATION  IN  BLACK  HILLS-1. 
Approximate  scene  dimensions  are  1166  ft  (355  m)  by  7212  ft 
(2198  m)  for  channel. 12  and  1166  ft  (355  m)  by  7146  ft  (2178  m) 
for  channels  2,  4,  5.  Each  sub-area  as  well  as  the  total  area 
have  been  histogrammed.  Histogram  plots  and  their  respective 
sub-areas  are  identified  with  the  following  key: 

□  Sub-area  1  +  Sub-area  4 

O  Sub-area  2  x  Sub-area  5 

A  Sub-area  3  O  Sub-area  6 
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AREA*  BLACK  HILLS 


00*00 
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It  BLACK  HILLS 


CALJB. PLATES 


AREA*  BLACK  HILLS  1 
LAMBDA*  8.0  TO  13.5 


BLACK  HILLS- 2 


Scene  Type  Forested  Mountains 

Date  of  Flight  22  July  1969 


Time 

of 

Flight 

1340  -  1342 

Altitude  (Ft) 

1500 

No. 

of 

Sub-Areas 

6 

No. 

of 

Data  Points 

468,915 

Channels 

5 

3 

7 

Wavelength  (um) 

1.0-1. 4 

1.5-1. 8 

2. 0-2. 6 

Resolution  (mr) 

In-Track 

6.6 

6 . 6 

6.6 

Cross-Track 

3.5 

3.5 

3.5 

Nadir  Pixel  Dimensions 
In-Track 

(m) 

3.017 

3.017 

3.017 

Cross-Track 

1.600 

1.600 

1.600 

Nadir  Ground  Sample 

Distance  (m) 

In-Track 

3.017 

3.017 

3.017 

Cross-Track 

1.143 

1.143 

1.143 

Line  Averaging  used  for  ALL  channels. 


*The  Black  Hills-2  data  were  collected  with  an  M-5  scanner  with 
calibration  reference  lamps  in  the  scanner  housing.  The  Black  Hills-1 
and  Black  Hills-2  data  were  collected  at  the  same  time  with  two  M-5 
scanners  in  the  aircraft. 
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MAP  OF  BLACK  HILLS -2 
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LINE  SCAN  IMAGES  PRODUCED  FROM  THE  VARIOUS 
INFRARED  CHANNELS  OF  BLACK  HILLS-2 


M 

r§|gg 

R 

gMflfr  tfk 'JB&f 

Pixel  1 
Pixel  161 


Line  2 


Line  364 


Line 


Pixel  484 

Pixel  645 
728 


SUB-AREAS  DEFINED  FOR  STATISTICS  GENERATION  IN  BLACK  HILLS-2  IMAGE. 
Uneven  areas  were  chosen  so  that  Areas  2  and  5  covered  the  ±20° 
range  suitable  for  correlation.  Approximate  scene  dimensions  are 
2418  ft  (737  m)  by  7206  ft  (2196  ra) .  Each  sub-area  as  well  as  the 
total  area  have  been  histogrammed.  Histogram  plots  and  their 
respective  sub-areas  are  identified  with  the  following  key: 

a  Sub-area  1  +  Sub-area  4 

O  Sub-area  2  x  Sub-area  5 

A  Sub-area  3  ❖  Sub-area  6 
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BLACK  HILLS- 2 
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AREA:  BLACK 


AAEAi  BLOCK  HILLS  2 

LAMBOft*  t.5~T0  1.8  A*  HERN  -  447.21 

TOTAL  AREA  ST. DEV.-  154.85 


mLIB.PlATfS 


I 


PISGAH  CRATER* 


Scene  Type 
Date  of  Flight 
Time  of  Flight 
Altitude  (Ft) 

No.  of  Sub-Areas 
No.  of  Data  Points 


Mountains 
30  October  1970 
0822  -  0827 
1000 
4 

650,076  for  channel  4 
637,632  for  channels  1,  2 


Channels 


4  12 


Wavelength  (pm) 


11.3-13.5  8.0-10.9  9.4-12.1 


Resolution  (mr) 
In-Track 
Cross-Track 


3.5  28  28 

3.5  21  21 


Nadir  Pixel  Dimension  (m) 

In-Track  1.067 

Cross-Track  1.067 


8.534  8.534 

6.401  6.401 


Nadir  Ground  Sample 
Distance  (m) 

In-Track  1.067  1.067 

Cross-Track  0.762  0.762 


1.067 

0.762 


No  Line  Averaging  was  used. 


*The  Pisgah  Crater  data  were  collected  with  an  M-5  scanner  with 
thermal  calibration  plates  in  part  of  the  scanner  f ield-of-view. 
The  11.3-13.5  urn  detector  and  the  8.0-10.9,  9.4-12.1  urn  detectors 
were  on  opposite  ends  of  the  scanner  and  are  not  in  spatial 
registration.  Hence,  spectral  correlation  coefficients  have 
not  been  determined  between  the  11.3-13.5  pm  data  and  either 
the  8.0-10.9  or  the  9.4-12.1  pm  data.  Histograms  and  correlation 
coefficients  for  the  8.0-10.9  and  9.4-12.1  pm  data  are  shown 
only  for  sub-areas  2,  3,  and  4.  Data  in  sub-area  1  were  lost  in 
the  analog-to-digital  conversion  process. 
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IMAGES  PRODUCED  FROM  THE  VARIOUS  INFRARED  CHANNELS  OF  PISGAH  CRATF.R 


PISGAH  CRATER 


SUB- AREA  1 


CHAMNt'l  S 


2.902«F*O2 


MEAN 


2.7(»fl8E  +  ''0 


TOTAL  PTS 


PISGAH  CRATER 
SUB-AREA  2 


ai  in; 
i 


l  .""f> 


.‘it.li.1 _ L  .£111! 


rK A 'i A  TI  ■; 
nr  a ‘‘ 

Ul  ,  JJF  V. 

1 1  *  T  A  |  MT. 


1  ...  ..  P 

n.f<o  i»r  m? _ j.gqon  ♦(?«» 
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4> 


PISGAH  CRATER 
SUB-AREA  3 

rnpwfcLM  InN _ i _ ? _ 

_ i _ 1 .  o  u  o _ 

_ 2 _ o.7Sb  i.onp 


CHANNEL  S 

1 

2 

u 

MEAN 

2.*e2tr+n? 

2>ftft5F.*02 

2.«*7'?r  +  n? 

SI.  DEV. 

2.07061*00 

?.P»  7ir  too 

1HTAI  PT3. 

_ 
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It  PISGflh  CRATER 
OA-  11.3  TO  13.5*41 


AREA i  PISGAH  CRATER 
LAMBOA-  11.3  TO  13.5 
CALIB. PLATES 
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AREA*  PISGAH  CRATER 
LAMBOA-  8.0  TO  10.9  nH 


RflEflt  PISGflH  CttflTEfl 
LflHBOfl®  9.4  TO  12.1#* 


AREA*  PISGAH  CRATtfl 
LAMBOA=  9.4  TO  12.1  pH 
CALIB. PLATES 


MONO  LAKE 


Scene  Type 
Date  of  Flight 
Time  of  Flight 
Altitude  (Ft) 

No.  of  Sub-Areas 
No.  of  Data  Points 


Mountains 
23  September  1968 
0952  -  0957 
4000 
4 

79,360  for  channels  2,  4,  5 
87,730  for  channel  20 


Channels 

2 

4 

5 

20 

Wavelength  (pm) 

1.0-1. 4 

2. 0-2. 6 

4. 5-5. 5 

8.0-13. 

Resolution  (mr) 

In- Track 

6.6 

6.6 

6.6 

6.6 

Cross-Track  3.5 

Nadir  Pixel  Dimension  (m) 

3.5 

3.5 

6.6 

In-Track 

8.046 

8.046 

8.046 

8.046 

Cross-Track 

4.267 

4.267 

4.267 

8.046 

Nadir  Ground  Sample 

Distance  (m) 

In- Track 

8.046 

8.046 

8.046 

8.046 

Cross-Track 

3.048 

3.048 

3.048 

3.048 

Line  Averaging  used 

for  channels 

2 ,  4 ,  and 

5  and  for 

channel  20 

independently. 


*The  Mono  Lake  data  were  collected  with  an  M-5  scanner  with 
thermal  calibration  plates  in  part  of  the  scanner  f ield-of-view. 
The  8.0-13.5  pm  detector  and  the  1.0-1. 4,  2. 0-2. 6,  4.5-5. 5  pm 
detectors  were  on  opposite  ends  of  the  scanner  and  are  not  in 
spatial  registration.  Hence,  spectral  correlation  coefficients 
were  not  determined  between  the  8.0-13.5  pm  data  and  the  1.0-1. 4 
2. 0-2. 6,  and  4. 5-5. 5  pm  data. 
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^t^SmSSeSSmS. 

IN  THE  MONO  LAKE  IMAGE.  Approximate  scene 
dimensions  are  7497  ft  (2285  m)  by  3110  ft 
(948  m)  for  channel  20  and  6758  ft  (2059  m) 
by  3110  ft  (948  m)  for  channels  2,  4,  5. 

Each  sub-area  as  well  as  the  total  area  have 
been  histogrammed .  Histogram  plots  and  their 
respective  sub-areas  are  identified  with  the 
following  key: 

□  Sub-area  1  +  Sub-area  4 

O  Sub-area  2  x  Sub-area  5 


A  Sub-area  3  O  Sub-area  6 


mW/CM2  STR  juM 


CALIBRATES 
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AREA*  MONO  LAKE 

LAMBOA-  2.0  TO  2.6  mH  MEAN  -  94.65 

TOTAL  AREA  ST. DEV.-  16.00 


RREflt  MONO  LflKE 


CAL IB. PLATES 


300.00  305 

KELVIN 
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CALIB.  PLATES 


WIENER  SPECTRA  AND  AREA/ INTENSITY 
STATISTICS  FOR  FLINT-1  AND  MILL  CREEK 


-  ■ 


JIXI  fPIXI  gOIXI  gOIXl  jOIXl  1  x-01 

(W/I)/2**ndWb)  A1ISN3G  d3M0d 
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SPATIAL  FREQUENCY  (1/M) 


cjOTXI  ipixi  £01X1  gOIXI  xOIXI  T  x-jOIXI  2jOIX1 

(W/I3  /Z**  ridWU)  JL1ISN30  y3M0d 


n-192 


flREfli  FLINT  1 
I NT RACK 


(W/I)/2**ndHU)  A1ISN3G  d3M0d 
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AREA*  MILL  CHEEK 
INTRACK 


323 
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FLINT  1 

RflOIRNCE  THRESHOLD  -  1.50CJ 
LflMBOfl-  1.0  TO  1.4  mH 
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FLINT-1 


DISTRIBUTION  OF  RfcCOOlZEO  PAOIANCE  THRESHOLDS  FOR  a  -  1.5 

BY  AREA 


SQUARE 

MFTERS 

FREQUENCY 

0.0 

TU 

100.0 

1012 

100.0 

TO 

2»0.0 

30 

200.0 

TO 

500. 0 

17 

500.0 

TO 

1000. 0 

4 

Inoo.o 

TO 

1500.0 

0 

1500.0 

TO 

2000.0 

0 

2ooo.o 

TO 

2500.0 

1 

2500.0 

TU 

3000.0 

0 

3ooo,o 

TO 

4000.0 

0 

flooo.o 

TU 

5oon. 0 

0 

5O00.0 

TU 

6000.0 

0 

6000.0 

TO 

0000. 0 

0 

0000 . 0 

TU 

1  0000.0 

0 

l  oooo.o 

TO 

15000.0 

0 

1 5(100.0 

TO 

200O0.0 

0 

20000.0 

TO 

'1 0  0  o  0 . 0 

0 

00000.0 

TO 

soooo. o 

0 

80000.0 

TO 

1601)00.0 

0 

OVER 

160000.0 

0 

TOTAL  NUMBER  OF  RADIANCE  TMRs  1064 
2340  FEATUHFS  «ITH  areas  LESS  Than  3. Si  HECTARES  were  ALSO  RE COGNIZED 


METERS 

BY  PERIMETER 

FFET 

FREQUENCY 

SHAPE 

BY  SHARE 

FACTUR  FRFQUE 

0 

TO 

50 

0 

Tf) 

164 

079 

0.0 

TO 

1.0 

2 

50 

TO 

100 

164 

TO 

320 

120 

1.0 

TO 

1.1 

0 

100 

TO 

150 

320 

TO 

44£ 

29 

1  .  1 

TU 

1.2 

43 

150 

TO 

200 

0  92 

Tn 

656 

17 

1 .2 

TO 

1.3 

41 

2(>0 

TO 

250 

656 

Tn 

820 

5 

1.3 

TO 

1.4 

94 

250 

TO 

300 

020 

TO 

904 

1 

1.4 

TP 

1.5 

97 

300 

TO 

350 

904 

TP 

1  140 

2 

J .  5 

TP 

1.6 

87 

350 

TO 

400 

1  140 

TO 

1312 

5 

1.6 

TP 

1.7 

«2 

400 

TO 

500 

1312 

TO 

1640 

2 

1.7 

TO 

1.8 

94 

500 

’0 

600 

1  640 

TO 

I960 

1 

1.0 

TO 

1.9 

53 

600 

TO 

700 

1900 

TO 

22  96 

2 

1.9 

TO 

2.0 

75 

700 

TO 

800 

2206 

TP 

2624 

1 

2.0 

TO 

2.2 

99 

800 

TO 

900 

2629 

TO 

2952 

0 

2.2 

TO 

2.4 

66 

900 

TO 

1000 

2952 

TO 

32«0 

0 

2.4 

Tn 

2.6 

52 

1000 

TO 

1200 

320  0 

TO 

3937 

0 

2.6 

TO 

2.0 

60 

1200 

TO 

1  900 

39  37 

Tf) 

0 

2.C 

TP 

3.0 

25 

1400 

TO 

1600 

4593 

TU 

5244 

0 

3.0 

TO 

3.5 

44 

1600 

TO 

2000 

5249 

TO 

6561 

0 

3.5 

TC. 

4.0 

22 

UVFR 

2000 

OvF.F. 

656  l 

0 

.  OVCR 

4.0 

20 
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FLINT-1 


* 

DISTRIBUTION  OF  PFCnGvlzEw  RADTANCE  THRESHOLDS  FOR  a  -  2.0 

MY  ARFA 


square 

METERS 

FREQUENCY 

0,0 

TO 

l 0  0 . 0 

167 

loo.o 

TO 

200.0 

2 

200.0 

TO 

5*10.0 

l 

500,0 

TO 

loon.o 

0 

1 ooo.o 

TO 

151*0.0 

0 

ISOo.o 

TO 

200  0.0 

1 

?ooo.o 

TO 

25<>0.0 

0 

2500.0 

TO 

3<>oo .  o 

0 

3o0o.o 

TU 

aiioo.o 

0 

0000. 0 

TO 

501*0.0 

0 

500o.o 

TU 

60 1*0.0 

0 

oooo.o 

TO 

*000.0 

0 

BoOo.o 

TO 

1  o  n  o  o .  o 

0 

10000. 0 

TO 

15900 .o 

0 

15000.0 

TO 

29000.0 

0 

20000.0 

TO 

9  0  (1  (*  0 , 9 

0 

OOOOtl.o 

TM 

09000.0 

0 

*0000. 0 

TO 

169090.0 

0 

OVER 

16901*0.0 

0 

TOTAL  NUMBER  OF  RADIANCE  Thr=  191 

759  FEATURES  wIT*  A»paS  LESS  Than  3.*l  HtCTARfcS  MERE  ALSO  RECOGNIZED 


METERS 

BY  PFRntTC* 

FFFT 

FRFQIJENCY 

SHAi»r 

9 Y  SHARE 

FACTOR  FRf.qOE 

0 

TO 

50 

0 

TO 

169 

176 

0.0 

TO 

1.0 

1 

50 

TO 

100 

169 

TO 

32* 

1  1 

i  .0 

TO 

t.l 

0 

190 

TO 

150 

326 

TO 

992 

2 

i.i 

TO 

!.? 

0 

150 

TO 

200 

902 

TO 

666 

1 

1.2 

TO 

1.3 

9 

200 

TO 

250 

656 

TP 

*20 

0 

1.3 

Trj 

1.0 

2* 

250 

TO 

300 

«2o 

TO 

9*4 

0 

1 .0 

TO 

1.5 

12 

300 

TO 

350 

9<ia 

TO 

1 1  T* 

0 

1.6 

TO 

1.6 

1* 

350 

TO 

990 

119" 

TO 

1312 

0 

1.6 

TP 

1.7 

17 

9  00 

TO 

500 

131? 

TO 

1 6  0  0 

0 

1.7 

TO 

1.* 

1* 

500 

TO 

600 

1699 

in 

196* 

0 

1.* 

TO 

1.9 

1 1 

600 

TO 

700 

196" 

TO 

2296 

0 

1.9 

TO 

2.0 

22 

700 

TO 

MOQ 

2296 

in 

26 20 

1 

2.0 

TU 

2.2 

20 

000 

TO 

900 

2629 

TO 

2*>52 

0 

2.2 

TO 

2.0 

10 

990 

TO 

1  090 

2952 

TO 

32*0 

0 

2.9 

TO 

2.9 

8 

1000 

TO 

1200 

32*0 

TO 

3937 

0 

2.6 

TO 

2.H 

3 

1200 

TO 

1900 

3937 

TO 

969  3 

0 

2.* 

TO 

3.0 

1 

19  00 

TO 

1600 

0593 

TO 

S2U9 

0 

3.0 

TO 

3.5 

2 

1600 

TO 

2000 

5209 

TO 

6‘ih  l 

0 

3.5 

TO 

9.0 

2 

OVER 

2000 

OVF.K 

656  1 

0 

OVCR 

9.0 

l 
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FLINT  I 

RADIANCE  THRESHOLD  -  2.50CT 
LAMBOA-  1.0  TO  1.4  wH 
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FLINT-1 


distribution  of  RF.cnr.iiizFn  radiance  turnouts  fox  a  -  2.5 

t'.V  area 

SQUARE  HETERS  FREQUENCY 


0,0 

TU 

too.o 

?o 

100.0 

TO 

200.0 

1 

200.0 

TO 

boo.O 

1 

50o. o 

TO 

t  0  0  o  .  0 

1 

I  0  0  0 . 0 

TO 

1500.0 

0 

1500.0 

TO 

2oiio.« 

0 

2(100,0 

TO 

2b<i0.o 

0 

2500.0 

TO 

3o»<>.  o 

0 

3<»0o.o 

TO 

0000.0 

0 

OoOo. 0 

TO 

'iOO !) . 0 

0 

5o0o.o 

TO 

6000.0 

0 

6000. 0 

TO 

0000,0 

0 

Hooo.o 

TO 

]  0  0  0  o ,  0 

0 

10O0O.O 

TO 

15000. o 

0 

JSoOo.o 

TP 

20000.0 

0 

20O0O.O 

TO 

0  o  0  0  o  .  0 

0 

noooo.o 

TO 

sono o.o 

0 

eoooo.o 

TO 

1 6  o  0  <’  0.0 

0 

nvF« 

160000,0 

0 

TUTAL  wij-BER 

OF  PaOIAMCE  TW«s 

32 

100  FCATURfS  a!  T  TH  AREAS  LESS  Than  3. HI  HECTARES  *tPF  ALSO  RECOGNIZE" 


BY  perimeter  BY  SHAPE 


metfrs 

feet 

FPEmiE'XY 

SHAPf 

FACTOR 

FREQUENCY 

0 

TO 

50 

0 

TO 

160 

20 

0.0 

TO 

1.0 

0 

50 

Tn 

100 

160 

TO 

320 

5 

1.0 

TP 

1.1 

0 

too 

TO 

150 

320 

rn 

0  92 

1 

1.1 

TO 

1.2 

1 

150 

TP 

200 

092 

TO 

656 

O 

1.2 

tn 

1.3 

2 

200 

Tn 

250 

656 

TO 

020 

1 

1.3 

TO 

1.0 

2 

250 

T  n 

300 

«2'» 

TO 

9«0 

0 

1.0 

Tn 

1.5 

0 

300 

TO 

350 

900 

TO 

1  1  OH 

1 

1  .5 

TO 

1.6 

2 

350 

Tn 

400 

1100 

TO 

1312 

0 

1  .6 

TP 

1.7 

1 

000 

to 

500 

1312 

TO 

1600 

0 

1.7 

TO 

l.» 

3 

SCO 

Tn 

600 

J  600 

TO 

1 96H 

0 

1.0 

Tn 

1.9 

t 

600 

Tn 

700 

i960 

TO 

2296 

0 

1  .9 

TO 

2.0 

a 

700 

TP 

000 

22*»(. 

Tn 

262o 

0 

2.0 

TO 

2.2 

2 

800 

TP 

900 

26?a 

TO 

2952 

0 

2.2 

TO 

2.0 

0 

900 

TP 

1000 

2952 

TP 

3200 

0 

2.0 

To 

2.6 

0 

1000 

TP 

1200 

3280 

TO 

3937 

0 

2.6 

Tr» 

2.0 

1 

1200 

Tn 

1000 

39  37 

TO 

0593 

0 

2.0 

Tn 

3.0 

1 

touo 

TP 

1600 

059  3 

TO 

52oo 

0 

3.0 

Tn 

3.5 

3 

1600 

Tn 

2000 

5209 

TO 

656  1 

0 

3.S 

TO 

0.0 

1 

OVER 

2000 

ovfp 

656  1 

0 

OvCR 

0.0 

0 

ru  ru  ty  ru 
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FLINT  1 

RflDIRNC E  THRESHOLD  -  3.000" 
LRH60R-  1.0  TO  1.4  m 
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FLINT- 1 


DISTRIBUTION  OF  RECOGNIZED  BtOT AHCE  THRESHOLDS  FOR  <j  =  3.0 

BY  AREA 

SQUARE  METERS  FREQUENCY 


O.o 

TO 

100.0 

16 

1  Oo. n 

TO 

2«0. 0 

0 

20o.o 

TO 

500. 0 

0 

QOo. n 

TO 

looo.o 

1 

luoo.o 

TO 

!5oo. 0 

0 

»50o. n 

TO 

2000.0 

0 

200o.o 

TO 

25»n  .0 

0 

250o. o 

TO 

30"0.0 

0 

3o0o,n 

TO 

4000. 0 

0 

4000.0 

TO 

5000.0 

0 

Sooo.o 

TO 

6000. 0 

0 

(tOUO.O 

TO 

flooo.o 

0 

8ooo.  n 

TO 

1 Pooo.o 

0 

10000. n 

TO 

15000.0 

0 

1 5  o  0  o .  o 

TO 

20000. 0 

0 

20ouo.o 

TO 

nopnn.o 

0 

4  00  00 . 0 

TO 

eooop.o 

0 

800  0  o.o 

TO 

1600O0.0 

0 

OVER 

160000.0 

0 

TOTAL  NUMBER 

(IF  PAOJANCt  THRs 

17 

00  FEATURFS  WIT”  AREAS  ir.«3  1man  3.81  HECURLS  *EPF  also  RECOGNIZED 


BY  PFRl“fcTtR  HY  SHAPE 


METERS 

teet 

FREQUENCY 

SHAPE 

FACTOR 

FREQUENCY 

0 

TO 

50 

ii 

TO 

164 

12 

0.0 

TO 

1.0 

0 

50 

TO 

100 

164 

T»J 

328 

4 

1.0 

TO 

1.1 

0 

100 

TO 

150 

328 

TO 

4  92 

0 

1.1 

TO 

1.2 

2 

150 

TO 

200 

49? 

in 

656 

0 

1  .2 

TO 

1.3 

1 

200 

TO 

250 

656 

TO 

820 

0 

1.3 

TO 

l.« 

3 

250 

TO 

300 

82" 

TO 

984 

0 

1.4 

Tn 

1.5 

1 

300 

TO 

350 

OAK 

TO 

1148 

0 

1.5 

To 

1.6 

3 

350 

TO 

400 

1148 

TO 

1312 

0 

1  .6 

TO 

1.7 

0 

400 

TO 

500 

1312 

TO 

1640 

1 

1.7 

TO 

1.8 

0 

500 

TO 

600 

1640 

TO 

I960 

0 

1.8 

10 

1.9 

0 

600 

TO 

700 

19*0 

To 

2296 

0 

1.9 

TO 

2.0 

1 

700 

TO 

800 

2296 

TO 

2624 

0 

2.0 

TO 

2.2 

2 

800 

TO 

900 

2624 

To 

2952 

0 

2.2 

TO 

2.“ 

2 

900 

TO 

1000 

2"52 

TO 

3289 

0 

2.4 

TO 

2.6 

0 

1000 

TO 

T  200 

320" 

TO 

3937 

0 

2.6 

TO 

2.8 

0 

1200 

TO 

1400 

3937 

TO 

4593 

0 

2.4 

TO 

3.0 

1 

1400 

TO 

1600 

459  3 

TO 

5249 

0 

3.0 

TO 

4.5 

0 

1600 

TO 

20«0 

52aO 

TO 

6561 

0 

3.5 

Tn 

4.0 

0 

OVER 

2000 

ovfr 

6561 

0 

OVER 

4.0 

1 
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FLINT  1 

TEMPERATURE  THRESHOLD  -  1.50C 
LAMBOA-  9.3  TO  11.7  mM 


11-204 


FLINT-1 


DISTRIBUTION  UF  RECOGNIZED  TfMP.  THRESHOLDS  FOR  a  »  1.5 


square 

HETFRS 

FREQUENCY 

0.0 

TO 

100.0 

1278 

tOO.n 

TO 

2<)0.0 

20 

200.0 

TO 

500. 0 

4 

500,0 

TO 

1000.0 

2 

looo.o 

TU 

15«»0.0 

0 

1500.0 

TO 

2000.0 

0 

2000.0 

TO 

25oo. 0 

0 

2500.0 

TO 

3ooo. 0 

0 

3uon.n 

TO 

floor). 0 

0 

QUO  0.0 

TO 

5000. 0 

0 

Sooo.o 

TO 

6000.0 

1 

boOO.O 

TO 

8000. 0 

0 

Sooo.o 

TU 

10000.0 

0 

10000.0 

TO 

15000.0 

0 

15000.0 

TU 

20000.0 

0 

20000.0 

TO 

0  o  0  n  0 , 0 

0 

4  0  0  0  0 . 0 

TO 

80000,0 

0 

80000.0 

TU 

lbOOOO.O 

0 

OVER 

1600 0 0.0 

0 

TOTAL  HUMBER 

UF  TEMP. 

T  hr  s  nos 

370  FEATURES  with  arfaS  LESS  THAN  3.81  hECIARES  WE»f  ALSO  RECOGnizEO 


BY  PERIMETER  HY  SHAPE 


METERS 

FFET 

FREQUENCY 

SHAPE 

.  FACTOR 

F«C JUfcNCY 

0 

TU 

50 

0 

TO 

164 

1254 

0.0 

TO 

1.0 

3 

50 

TO 

too 

164 

TO 

328 

42 

J  .0 

TU 

1.1 

3 

100 

TU 

ISO 

320 

TU 

492  * 

5 

1.1 

TO 

1.2 

104 

ISO 

TO 

200 

492 

TU 

65r> 

2 

1  .2 

TO 

1.3 

186 

20  0 

TU 

2S0 

656 

TU 

820 

0 

1.3 

TP 

1.4 

S  39 

2S0 

TU 

300 

820 

TO 

98fl 

1 

1.4 

TO 

1.5 

238 

300 

TU 

350 

98U 

TO 

1148 

0 

1.5 

TO 

1.6 

167 

3S0 

TU 

900 

1  IflO 

TP 

1312 

0 

1.6 

TP 

1.7 

96 

400 

TU 

500 

1312 

TO 

1640 

0 

1.7 

TO 

l.« 

56 

500 

TU 

600 

160  0 

TU 

1  068  ‘ 

0 

1.8 

TO 

1.9 

31 

600 

TU 

700 

1«>60 

TO 

2296 

0 

1  .9 

TO 

2.0 

26 

700 

TU 

800 

2296 

TO 

2624 

0 

2.0 

TU 

2.2 

33 

800 

TU 

900 

262  o 

TU 

2952 

l 

2.2 

TO 

2.4 

lu 

900 

TU 

1000 

2952 

TU 

3280 

0 

2.4 

TO 

2.6 

4 

1000 

TU 

1200 

3280 

TO 

39.37 

0 

2.6 

TO 

2.6 

3 

1200 

TU 

1  flOO 

3937 

TO 

4593 

0 

2.8 

TO 

3.0 

4 

lflOO 

TU 

1600 

0593 

TO 

5249 

0 

3.0 

TO 

3.5 

l 

1600 

TU 

20  00 

5209 

TU 

6-561 

0 

3.5 

TO 

4.0 

1 

OVER 

2000 

UVFR 

6561 

0 

OVER 

4.0 

0 

11-205 
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FLINT  1 

TEMPERATURE  THRESHOLD  -  2.00CT 
LflMBOR-  9.3  TO  11.7  wM 


11-206 


FLINT-1 


BY  AREA 


SQUARE  METFRS  FREQUENCY 


0.0 

TU 

100.0 

870 

100.0 

U) 

2»0.0 

8 

200. 0 

TU 

500.0 

a 

500. o 

TO 

loo o.o 

2 

1000.0 

TO 

15»o. 0 

0 

1500.0 

TO 

?0oo. 0 

0 

2O0O.0 

TU 

2500.0 

0 

2500.0 

TO 

3ooo.O 

0 

3000.0 

TU 

0000.0 

0 

4000.0 

TO 

•1000.0 

1 

5000.0 

TO 

60  0  0,0 

0 

6000.0 

TO 

8000 . 0 

0 

8000,0 

TO 

1  0  0  0  0 . 0 

3 

1  0  0  0  0 . 0 

TO 

15000.0 

0 

15000.0 

TO 

20n('0.0 

0 

20000.0 

TO 

n  o  n  on .  o 

0 

ooooo.o 

TO 

80  0  0  0.0 

0 

aouoo.o 

TO 

1600(10.0 

0 

nvEH 

160000.0 

0 

TOTAL  NUMBER  UF  TE^P.  TMRs  88< J 


293  FEATURES  wItm  ARFAS  LESS  Than  3.61  HECTARES  WERE  ALSO  RECOGNIZED 


BY  Ph  R I ME  TLR 

METERS  FFET  FREOUENCY 


0 

TO 

50 

(1 

TO 

16« 

870 

50 

rn 

100 

164 

TP 

328 

10 

100 

TO 

150 

326 

To 

492 

2 

150 

TO 

200 

492 

TP 

656 

t 

200 

TO 

250 

656 

ip 

820 

l 

250 

TO 

300 

820 

Tn 

984 

0 

300 

TO 

350 

984 

10 

1  1  48 

0 

350 

TO 

UOO 

1108 

Tn 

1312 

0 

UOO 

TO 

500 

1312 

TO 

1  o'JO 

0 

5oO 

TO 

600 

1 60  n 

Tn 

1908 

0 

600 

TO 

700 

J9hfl 

TO 

2296 

0 

700 

TO 

HOC 

2296 

TO 

26?J 

0 

800 

TO 

900 

2624 

TO 

2952 

0 

9uO 

TD 

1000 

2952 

TP 

3260 

0 

1000 

TO 

1200 

3280 

to 

3937 

l 

1200 

TO 

1400 

3R37 

TO 

4593 

0 

1000 

TO 

1600 

059  3 

T'l 

5249 

0 

1600 

TO 

2000 

5209 

T'l 

656  1 

0 

OVF.R 

2000 

Over 

6561 

0 

rtV  SHAPE 


SHAPE.  FACTOR  FRfjUENCY 


0.8 

Tn 

1.0 

4 

J  ,  8 

TP 

1.1 

1 

1.1 

TO 

1.2 

63 

1.2 

TO 

1.3 

133 

1.3 

TO 

1.4 

172 

1.0 

TP 

1.5 

104 

1  .5 

T !  1 

1.6 

116 

1.6 

TO 

1.7 

93 

1.7 

TO 

1.8 

07 

1.8 

TP 

1.9 

42 

1.9 

TO 

2.0 

37 

2.8 

TO 

2.2 

20 

2.2 

T*i 

2.0 

10 

2.4 

Tn 

2.6 

3 

2.6 

TP 

2.6 

2 

2.8 

To 

3.0 

1 

3.0 

TP 

3.5 

0 

3.5 

TO 

4.0 

0 

pvcr 

4.0 

1 

DISTRIBUTION  OF  RECOGMZEO  TEMP. 


thresholds  FOR  a  »  2.0 


11-207 


FLINT-1 


DISTRIBUTION  OF  RFCl)GtJl7EO  TEMP.  THRESHOLDS  FOR  o  -  2.5 

nv  ARE* 


SQUARE 

;  METERS 

FREQUENCY 

0.0 

TO 

l«0.0 

271 

100. « 

TO 

2')  0.0 

4 

200.0 

TO 

5«»o  .  0 

a 

500.0 

TO 

looo.o 

1 

1600.0 

TO 

1500.0 

0 

1500.0 

10 

2000.0 

0 

2000.0 

TO 

?5'>o.O 

0 

2500.0 

TO 

3non.o 

0 

3000.0 

TO 

0  000. 0 

1 

4u0n.0 

TU 

8000. 0 

0 

50 00.0 

TO 

6000.3 

0 

6000.0 

TO 

flooo.o 

0 

8000.0 

TO 

10000. 0 

0 

1  0000.0 

TO 

1 500  0 . 0 

0 

15000.0 

TO 

20000.0 

0 

20000.  o 

TO 

onnoo.o 

0 

00000.0 

TO 

80000.0 

0 

80000.0 

TO 

16000 0.0 

0 

nvCR 

16000.0 

0 

TOTAL  NUM8FR  OF  TF«P.  TH«s  281 

202  FEATURES  *ITM  AREAS  LESS  THAN  3.81  HECTARES  *E  RE  ALSO  RECOGNIZED 


BY  PERIMETER  8 Y  SHAPE 


MCTERS 

FFET 

FREQUENCY 

SHAPE 

FACTUR 

FRt  JIJENCY 

0 

TO 

50 

0 

TO 

164 

269 

0.0 

TO 

t.o 

0 

50 

TO 

100 

1  6  M 

Tn 

328 

5 

1.0 

TO 

1.1 

1 

100 

TO 

150 

328 

TO 

4  9,? 

4 

1.1 

TO 

1.2 

27 

150 

TO 

200 

M92 

TO 

656 

1 

1.2 

TO 

1.3 

33 

200 

TO 

250 

656 

TO 

620 

l 

1.3 

TO 

1.4 

44 

250 

TO 

300 

820 

TO 

994 

0 

l.u 

Tn 

1.5 

29 

300 

TO 

350 

984 

TO 

1  1  48 

0 

1.5 

TO 

1.6 

51 

350 

TO 

MOO 

1  l'l8 

To 

1312 

0 

1.6 

TO 

1.7 

27 

400 

TO 

500 

1312 

TO 

1640 

0 

1 .7 

TO 

1.8 

26 

500 

TO 

600 

16M0 

TO 

1968 

0 

1." 

TO 

1.9 

11 

600 

TO 

700 

i960 

TO 

22=>6 

0 

1.9 

TO 

2.0 

1 1 

700 

TO 

800 

2290 

TO 

2024 

0 

2.0 

TO 

2.2 

12 

800 

TO 

900 

2624 

Tn 

2952 

1 

2.2 

TO 

2.4 

4 

900 

TO 

1  00U 

295? 

TO 

32«0 

0 

?.'» 

TO 

2.6 

0 

1000 

TO 

1200 

3280 

TO 

3937 

0 

2.6 

Til 

2.0 

1 

1200 

TO 

t  /■  0  0 

3937 

TO 

4593 

0 

■2.8 

TO 

3.0 

1 

I'lOO 

To 

1600 

45R3 

TO 

5249 

0 

3.0 

TU 

3.5 

0 

1600 

TO 

2000 

5?4  9 

TO 

656  1 

0 

3.5 

TO 

4.0 

0 

OVER 

2000 

over 

6561 

0 

OVER 

4.0 

1 

645 


FLINT  1 

TEMPERATURE  THRESHOLD  »  3.00(3" 
LAMBDA-  9.3  TO  11.7  mM 


11-210 


FLINT-1 


DISTRIBUTION  OF  REC'lDclZFO  TE-P.  THRESHOLDS  FOR  a  =  3.0 

BY  *PfA 

SQUARE  METFRS  FREOuE'.CY 


o.n 

TO 

100. 0 

48 

100.0 

TO 

200.0 

2 

200.0 

TU 

50(1.0 

1 

500.0 

TO 

1  noo.  r 

1 

l  000.0 

TO 

1500.0 

0 

I50O.0 

TU 

20"0.0 

0 

2o0o.n 

TO 

25»o.o 

0 

2500.0 

TO 

3000. 0 

1 

3000.0 

TU 

0  o  o  o  ,  o 

0 

4  o  OO.o 

TO 

'">000.0 

0 

5oco.o 

TO 

6000. 0 

0 

bo 00.0 

TO 

flood. 0 

0 

BbOo.o 

TU 

1  o  0  o  o  .  0 

0 

lbooo.o 

TO 

1  r>  0  o  0 . 0 

0 

15000.0 

TO 

20000.0 

0 

20O00.0 

TU 

n  o  o  o  o .  0 

0 

ooooo.n 

TO 

B00O0.0 

0 

SOUOO.O 

TO 

160000.0 

0 

OVER 

160000.0 

0 

TOtAL  NUMBER  OF  TFhP.  THRs  S3 

49  FEATURES  /i  I T  h  ARFaS  LESS  Thau  S.a|  HECTaRFS  "ERE  ALSO  REC'IGhIZEO 


METFRS 

BY  PERIm 

FEET 

E  T  C  R 

FREQUENCY 

SHAPE 

BY  S 

FAC1UR 

0 

TO 

50 

TO 

164 

Uo 

0.0 

TO 

1.0 

50 

TO 

mo 

164 

TU 

32fl 

3 

1.0 

TO 

1.1 

loo 

TO 

150 

32A 

TO 

4  92 

0 

1.1 

TO 

1.2 

150 

TO 

200 

402 

TO 

656 

3 

1.2 

TO 

1.3 

20  0 

TO 

250 

656 

TO 

820 

0 

1.3 

TU 

1.4 

250 

TO 

300 

fl.?o 

T(1 

904 

0 

1  .4 

TO 

1.5 

300 

TO 

350 

9rt'J 

TO 

1  1  16 

0 

1.5 

TO 

1.6 

350 

TO 

400 

1  14" 

to 

1312 

0 

1.6 

TO 

1.7 

000 

10 

500 

1  31? 

TU 

1640 

0 

1  .7 

T.J 

1.8 

500 

TO 

600 

164" 

Tn 

196H 

0 

1  .8 

TO 

1.9 

600 

TO 

700 

I960 

TO 

22=>6 

0 

1.0 

TO 

2.0 

700 

TO 

A00 

2?96 

TO 

2624 

0 

2.0 

TO 

2.2 

800 

TO 

900 

2l.?4 

TO 

2952 

1 

2.2 

TO 

2.4 

900 

TO 

moo 

2952 

TO 

32*  > 

0 

2.4 

TO 

2.6 

1000 

TO 

1200 

321" 

Tn 

39  37 

0 

2.6 

TO 

2.8 

1200 

TO 

140(1 

3<»37 

TO 

459  3 

0 

2.  A 

TO 

3.0 

14  00 

TO 

toOO 

'ISO  3 

Tn 

52  19 

0 

3.0 

TO 

3.5 

1600 

TO 

2000 

5  2  '*  9 

TO 

656  l 

0 

3.5 

TU 

4.0 

OVER 

2000 

OvE  K 

656 1 

0 

OVC  R 

4.0 

FRE.1UEnCY 

0 
I 
3 
1 1 


»  -O  l/\  «  <\l  IT  O 


MILL  CREEK 


DISTRIBUTION  OF  RECUUMZED  RADIANCE  THRESHOLDS  FOR  CT  -  1.5 

BY  AREA 

SQUARE  METERS  FREQUENCY 


0.0 

TO 

100. 0 

232 

100.0 

1U 

200. 0 

05 

200.0 

TU 

500.6 

25 

500.0 

TO 

I  000.0 

10 

1000. 0 

TO 

1SOO.O 

5 

1500.0 

TU 

?  000.0 

3 

PdOO.O 

TO 

2600.0 

0 

P 500. 0 

TO 

3000,0 

0 

3000.0 

TO 

0  0  0  0 . 0 

2 

oooo.o 

TO 

50  0  0.0 

3 

5000.0 

TO 

6000.0 

1 

6000.0 

TO 

0000. 0 

2 

s  0  0  0 . 0 

TO 

10000.0 

1 

1 0000.0 

TO 

1  5  0  0  0 , 0 

1 

1 5  0  0  0  •  0 

TO 

20000.0 

0 

20O0O.0 

TU 

00000. 0 

2 

OOoOO.O 

TO 

oonno.o 

1 

SOoOO.O 

TU 

160000. 0 

0 

OVER 

1600(10. 0 

0 

TOTAL  NUMBER  OF  RADIANCE  THRs  3:47 


BY  perimeter  BY  SHAPE 

METERS  FEET  FREQUENCY  SHAPE  FACTOR  FREQUENCY 


0 

TU 

50 

0 

TP 

160 

229 

o.n 

TO 

1.0 

6 

50 

TU 

100 

160 

TP 

328 

58 

1.0 

TO 

1.1 

2 

100 

TU 

150 

328 

TO 

09  2 

16 

1.1 

TO 

1.2 

105 

150 

TU 

200 

092 

TO 

656 

6 

1.2 

TO 

1.3 

27 

200 

TU 

250 

656 

TP 

620 

7 

1.3 

TO 

l.o 

03 

250 

TU 

300 

8  20 

TU 

980 

3 

l.o 

TO 

1.5 

10 

300 

TU 

350 

9«fl 

TU 

1108 

2 

I.S 

TU 

1.6 

1  1 

350 

TU 

0  00 

1108 

TO 

1312 

1 

1.6 

TU 

1.7 

22 

000 

TU 

500 

131? 

TP 

1600 

2 

1.7 

TO 

1.8 

15 

500 

TO 

600 

160  0 

TU 

1968 

0 

1.8 

T(1 

1.9 

7 

600 

TU 

700 

19*8 

TP 

2296 

2 

1.9 

TO 

2.0 

12 

700 

TU 

800 

2?  96 

TO 

26?o 

0 

2.0 

Tn 

2.2 

9 

800 

TO 

900 

262« 

TP 

2952 

0 

2.2 

TU 

2.0 

3 

900 

TO 

tooo 

295? 

Tn 

3280 

0 

2.0 

TU 

2.6 

b 

1000 

TU 

1200 

3280 

TP 

3937 

0 

2.6 

TO 

2.8 

2 

1200 

TU 

1000 

3937 

TP 

059  3 

0 

2.8 

TO 

3.0 

2 

1000 

10 

1600 

0593 

TP 

52«9 

2 

3.0 

TO 

3.5 

1 

1600 

TU 

2000 

52oo 

TO 

6561 

0 

3.5 

TU 

<1.0 

3 

OVER 

2000 

OvFH 

656  1 

1 

OVER 

O.o 

7 

11-213 


MILL  CREEK 


o 

DISTRIBUTION  of  PECoGn12ED  padIA^CF  THRESHOLDS  FOR  O  -  2.0 

BY  ARFA 


SQUARE 

METFUS 

FREQUENCY 

0,0 

TO 

100. 0 

165 

100.0 

TO 

200.0 

20 

200.0 

TO 

5oo.O 

26 

500.0 

TO 

1000.0 

7 

1000.0 

TO 

1500.0 

3 

1500,0 

TO 

20  0  0.0 

2 

2o0o,n 

TO 

25  o  o.o 

1 

2500.0 

TO 

3oi)o.O 

0 

3000. 0 

TO 

0000.0 

1 

0 1)00.0 

TO 

5non. 0 

0 

5ooo.o 

TO 

f>  0  0  o  •  o 

l 

(tbOO.O 

TO 

AlHIO.O 

0 

Sooo.o 

TO 

1  on  oo.O 

0 

1  0O0O.0 

TO 

15ono,o 

0 

lSuoo.o 

TO 

20000.0 

1 

20000.0 

TO 

Aorioo.n 

2 

OOoOO.O 

TO 

aonoo.o 

0 

B  0  0  0  0 , 0 

TO 

160000.0 

0 

OVER 

160000,0 

0 

TOTAL  NUMBER  OF  p  *ni AvCE  Thp=  229 


BY  PEPI MtTER 

9 Y  SHAPF 

HCTFR3 

FEET 

FREQUENCY 

SHAPF 

FACTOR 

FREQUEI 

0 

TT 

50 

0 

TO 

169 

156 

0.0 

TCI 

1.0 

1 

50 

TO 

100 

160 

Tfl 

32h 

39 

1.0 

TO 

1.1 

0 

too 

TO 

150 

320 

TO 

992 

15 

1.1 

TO 

1.* 

93 

150 

Tn 

200 

«92 

Tn 

656 

7 

1.2 

TO 

1.3 

13 

200 

TO 

250 

656 

Tn 

620 

3 

1.3 

TO 

l.« 

26 

250 

TO 

300 

02n 

TO 

909 

2 

1.9 

Tn 

1.5 

16 

300 

Tn 

350 

9HU 

Tn 

1196 

1 

1.5 

Til 

1.6 

15 

350 

TO 

«O0 

ii*» 

TO 

1  312 

0 

1.6 

Tfl 

1.7 

23 

000 

TO 

500 

1312 

TO 

1690 

1 

1.7 

TO 

1.8 

5 

500 

TO 

600 

16 '10 

TO 

1966 

2 

1.9 

in 

1.9 

4 

600 

TO 

700 

196" 

Tfl 

2296 

1 

1.9 

TO 

2.0 

4 

700 

TO 

600 

2296 

TO 

2629 

1 

2.0 

TO 

2.2 

10 

BOO 

TO 

900 

?62'» 

TO 

2952 

0 

2.2 

Tn 

2.9 

4 

900 

TO 

1000 

2952 

Tfl 

3260 

0 

2.0 

TO 

2.6 

3 

1000 

TO 

1200 

3260 

TO 

3937 

0 

2.6 

TO 

2.8 

3 

1200 

TO 

10  00 

393  7 

T.l 

9593 

0 

2.9 

Tn 

3.0 

0 

1000 

TO 

1600 

059  3 

TO 

529  9 

1 

3.0 

TO 

3.5 

3 

1600 

TO 

2000 

520  9 

Tfl 

6561 

1 

3.5 

TO 

9.0 

3 

OVER 

2000 

nvFH 

6SM 

2 

OVER 

9.0 

3 

11-215 


MILLCREEK 


RflOIflNCE  THRESHOLD  -  2.50<J* 
LflHBOR-  1.0  TO  1.1 


HILL  CREEK 


distribution  of  recosnizfo  raotance  thresholds  for  a  -  2.5 

RV  ASF* 


SQUARE 

;  PETERS 

FREQUENCY 

0.0 

TO 

100. 0 

138 

100.0 

TO 

200.0 

21 

200.0 

TO 

500. o 

15 

500.0 

TU 

1000.0 

2 

1  000.0 

TO 

1500. 0 

0 

1500.0 

TO 

2000,0 

a 

2000.0 

TO 

2500.0 

o  • 

2500.0 

TU 

3000.0 

0 

3000.0 

TO 

a  do  o.O 

0 

9ooo.  o 

TO 

5000.0 

l 

5000.0 

TO 

6  0  0  o  .  0 

0 

61)00.0 

TO 

8000. 0 

l 

8000. 0 

TU 

1  0000. 0 

l 

1 Oooo.n 

TO 

15000.0 

0 

15O00.0 

TU 

20000. 0 

0 

20000.0 

TO 

00000.0 

1 

noooo.o 

TO 

80000.0 

0 

80000. 0 

TO 

160000. 0 

0 

OVER 

160  oil  0.0 

0 

TOTAL  NUMBER  OF  RADIANCE  THR  =  180 


BY  PERI*t»CR  BY  SHAPE 


KFTFRS 

FEET 

FREQUENCY 

SHAPE 

FACTOR 

FRE  JlltNCY 

0 

TO 

50 

0 

TU 

169 

139 

0.0 

TO 

1.0 

0 

50 

in 

100 

169 

TM 

328 

26 

1.0 

TO 

1.1 

0 

too 

Tn 

150 

328 

Til 

992 

6 

1.1 

TO 

1.2 

86 

150 

TO 

200 

992 

TO 

656 

6 

1  .2 

TO 

1.3 

3 

200 

TO 

250 

656 

Til 

820 

1 

1.3 

TO 

1.4 

20 

250 

TO 

300 

820 

Tn 

989 

0 

1  .a 

TO 

1.5 

9 

300 

TO 

350 

989 

TO 

1198 

3 

1.5 

TO 

1.6 

7 

350 

TO 

ooo 

1198 

TO 

1312 

0 

1.6 

Tn 

1.7 

1 1 

000 

TO 

500 

1312 

Til 

169  9 

0 

1.7 

Tn 

1.3 

10 

500 

TO 

600 

169  0 

TO 

1968 

0 

1.8 

Tn 

1.9 

9 

60  0 

TO 

700 

1  RhB 

TO 

229fl 

0 

1.9 

Tn 

2.0 

9 

700 

TO 

600 

2296 

TO 

26?  9 

0 

2.0 

TO 

2.2 

6 

800 

TO 

900 

2629 

TO 

2952 

0 

2.2 

TO 

2.9 

6 

ROO 

TO 

1000 

2 952 

TU 

3281) 

1 

2.9 

TO 

2.6 

0 

1000 

TO 

1200 

3280 

TO 

3937 

0 

2.6 

TO 

2.8 

1 

1200 

TO 

1  900 

3937 

Til 

959  i 

0 

2.8 

TO 

3.0 

l 

tooo 

TO 

1600 

9595 

TO 

5299 

2 

3.0 

TO 

3.5 

2 

1600 

TO 

2000 

5299 

TU 

6561 

0 

3.5 

TO 

9.0 

2 

OVER 

2000 

OVFW 

65h  l 

1 

nvfs 

9.0 

3 

11-217 


MILLCREEK 

RADIANCE  THRESHOLD  -  3.000" 
LAMBDA-  1.0  TO  1.4  mH 


11-218 


MILL  CREEK 


* 

DISTRIBUTION  op  PEC0GH1ZEO  RADIANCE  THRESHOLDS  FOR  a  *  3.0 

HY  AREA 

SQUARE  METFRS  FREQUENCY 


0.0 

TO 

100.0 

173 

100.0 

TO 

200.0 

7 

200.0 

TO 

500. 0 

7 

50o.  o 

Til 

1000. 0 

3 

1 000.0 

TO 

1500.0 

0 

1500.0 

TO 

2000.0 

0 

2000.0 

TO 

250i).  0 

0 

2500.0 

TO 

3000. 0 

0 

3000.0 

TO 

4000.0 

0 

BoOO.O 

TO 

5000. o 

0 

5ooo.o 

TO 

6000,0 

0 

60  00.0 

TO 

hooo.o 

0 

8000.  o 

TO 

1  8000. 0 

0 

1  OUOO.O 

TU 

1 5000.0 

0 

15000.0 

TO 

20000,0 

0 

20000.0 

TO 

00000.0 

0 

00000.0 

TO 

90000.0 

0 

80000.0 

Til 

1601)  (10, 0 

0 

OVER 

1 6  o  o  o  0 , 0 

0 

TOTAL  NUMBER  OF  RAOIAnCE  THR  =  190 


BY  perimeter  BY  SHAPE 


METERS 

FFF.T 

FREQUENCY 

SHAPE 

FACTOR 

FREQUENCY 

0 

TH 

50 

0 

TO 

164 

169 

0.0 

TO 

1.0 

0 

60 

TO 

100 

164 

TO 

328 

1  1 

1.0 

TO 

1  .  1 

0 

100 

TO 

150 

328 

TO 

«4J 

4 

1.1 

TO 

1.2 

128 

150 

Tn 

2  00 

492 

TO 

656 

5 

1.2 

TO 

1.3 

3 

200 

TH 

250 

656 

10 

820 

0 

1.3 

TO 

1.4 

14 

260 

Tn 

300 

820 

TO 

984 

1 

1.4 

TO 

1.5 

2 

300 

m 

360 

984 

TO 

1148 

0 

I  .5 

TO 

1.6 

14 

360 

TO 

400 

1148 

TO 

1312 

0 

1  .6 

TO 

1.7 

14 

000 

TO 

500 

1312 

TO 

1640 

0 

1.7 

TO 

t.« 

4 

500 

TO 

600 

1640 

TO 

1968 

0 

1.0 

TO 

1.9 

2 

6  00 

TO 

700 

t  9b  A 

TO 

22°6 

0 

1.9 

TO 

2.0 

1 

700 

TO 

800 

2296 

TO 

2624 

0 

2.0 

TO 

2.2 

2 

BOO 

TO 

900 

2624 

TO 

2952 

0 

2.2 

TO 

2.4 

2 

900 

TO 

1000 

2952 

TO 

32"0 

0 

2.4 

TO 

2.6 

0 

1000 

TO 

1200 

328'* 

TO 

3937 

0 

2.6 

TlJ 

2.0 

2 

1200 

TO 

1000 

39  37 

TO 

4593 

0 

2.0 

Tn 

3.0 

0 

10  0  0 

TO 

1600 

4593 

TO 

520  9 

0 

3.0 

TO 

3.5 

2 

1600 

TO 

2000 

5244 

Til 

65b  l 

0 

3.5 

TO 

4.0 

0 

OVER 

2000 

OVER 

o5b  l 

0 

over 

4.0 

0 

( 


11-219 


MILL  CREEK 

TEMPERATURE  THRESHOLD  =  1. 

LflMBOfi*  9.3  TO  11.7** 


MILL  CREEK 


DISPOSITION  t|F  F'ECOGM  ZED  TE^P.  T‘«PC SH'ilos  FOR  a  =  1.5 

R  Y  AREA 


SQUARE 

.  hFTfrs 

FRFDIIENCY 

0.0 

10 

100.0 

568 

100. 0 

TO 

2oo. 0 

70 

coo  .o 

TO 

500.0 

5a 

500.0 

10 

1 000.0 

21 

luon.n 

Tu 

1500.0 

n 

1 5oo.  n 

10 

2  o  o  o .  r 

3 

CoOo.o 

TO 

25oo.o 

I  . 

2500.0 

TU 

3000.0 

0 

5o0o.o 

TO 

0  (1  0  0 . 0 

o 

o  o  0  o .  o 

TO 

■iooo.o 

0 

5oo o.O 

TO 

6000.0 

l 

6000.0 

TO 

8ooo.  0 

0 

8000.0 

1.J 

1  0  1!  M  I)  .  f. 

1 

10O0O.O 

hi 

15000.0 

4 

l5ooo.o 

10 

20000.0 

0 

20000. 0 

TO 

'IPOOp.O 

0 

OOOOO.O 

10 

80000.0 

1 

BOoOo.O 

Tu 

160000.0 

0 

nvFR 

160  (iOn.O 

0 

TOTAL  HgMHER 

UP  TE“P. 

T'8K=  736 

eY  PERr'f.TFR  BY  SHAPE 


hftfr 

S 

FFET 

FREQUENCY 

SHAPE 

FACTOR 

FkEnUEUCY 

0 

TO 

50 

0 

TCi 

164 

553 

0.0 

TO 

1.0 

2 

50 

in 

100 

164 

Tn 

328 

91 

1.0 

TO 

l  .1 

0 

IDO 

TO 

150 

32" 

TO 

aoj 

37 

1.1 

Tn 

1.2 

378 

150 

TO 

200 

092 

TO 

656 

20 

1.2 

TO 

1.3 

28 

?II0 

n 

250 

656 

TO 

820 

10 

1.3 

TO 

1.9 

91 

250 

TO 

300 

820 

Tf) 

9Aa 

3 

1.4 

TO 

1.5 

29 

300 

to 

350 

980 

TO 

1  1 

4 

1.5 

TO 

1.6 

25 

350 

TO 

«00 

114" 

T(’ 

1312 

3 

1.5 

TO 

1.7 

37 

400 

TO 

500 

1312 

TO 

1640 

3 

1.7 

TO 

1.0 

31 

500 

To 

600 

16'iu 

T(1 

1968 

1 

1." 

Tn 

1.9 

24 

6o  0 

TO 

740 

1  9(,0 

Tf. 

2296 

2 

1.9 

Tn 

2.0 

15 

700 

To 

800 

2296 

TO 

2624 

3 

2.o 

Tn 

2.2 

28 

OoO 

TO 

900 

26  2  « 

Tf! 

2952 

l 

2.2 

Tn 

2.4 

17 

900 

To 

1000 

2°52 

Ifi 

3280 

0 

?.« 

Tf) 

2.6 

10 

1000 

TO 

1240 

3  2  on 

TO 

3937 

0 

2.6 

Tn 

2.8 

2 

1200 

To 

1  400  ' 

39  37 

TO 

4593 

0 

2.« 

TO 

3.0 

4 

1400 

TO 

1  600 

4593 

TO 

5?49 

0 

3.n 

10 

3.5 

6 

1600 

TO 

2000 

5249 

Tf. 

656  l 

0 

3.5 

TO 

4.0 

4 

OVER 

2000 

OvF  R 

65fcl 

5 

OVER 

4.0 

5 

AO-A07?  5e*  ENVIRONMENTAL  RESEARCH  INST  OF  MICHIGAN  ANN  ARBOR  IN— ETC  F/G  17/5 
STATISTICAL  ANALYSIS  OF  TERRAIN  BACKGROUND  MEASUREMENTS  DATA* (U) 

MAR  77  R  SPELLICY  »  J  BEARD  t  J  R  MAXWELL  N00123-76-C-0708 
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ADA 
077  584 


646 


MILLCREEK 

TEMPERATURE  THRESHOLD  -  2.00(3" 
LAMBDA-  9.3  TO  11.7  «M 


11-222 


MILL  CREEK 


i 


DISTRIBUTION  OF  PECOGuI  ZEI.'  TO'P.  THRESHOLDS  FOR  a  -  2.0 


HY  AkEA 


SOUARE  METERS 


FPEP'iENCV 


0,0 

TO 

t'10.0 

383 

100,0 

TO 

200.0 

37 

200,0 

TO 

5l'0.0 

28 

500.0 

TO 

1000.0 

10 

1000. 0 

TO 

15»0.0 

3 

1500.0 

TO 

20110.0 

2 

?uoo.n 

TO 

2500.0 

2 

250  0. 0 

TO 

3oi  >0.0 

1 

3000,0 

TO 

U0O0.0 

1 

oooo.o 

TO 

5000.0 

2 

5000.0 

TO 

ooon.O 

0 

6000.0 

TO 

8  0 «'  0 . 0 

1 

8000. 0 

TO 

1 nooo.o 

0 

toooo.n 

TO 

15oon.o 

0 

15000.0 

TO 

20000.C 

0 

20000.0 

TU 

00000. n 

1 

40000." 

TO 

80000. 0 

0 

80000.0 

Til 

160 000.0 

0 

OVER 

164000.0 

0 

TOT Al.  MUMMER 

OF  TF“F. 

THPs 

#7  1 

BY  PEPIf*ETEP 


HY  SHAPE 


METER 

5 

FEET 

FREQUENCY 

SHAPE 

FACTUR 

FPEUOfcUCY 

0 

TO 

50 

0 

TO 

164 

373 

0.0 

TO 

1.0 

0 

50 

Tn 

too 

1  60 

TO 

328 

49 

1.0 

TO 

1.1 

0 

too 

TO 

150 

324 

TO 

492 

20 

1.1 

Tn 

1.2 

255 

150 

TO 

200 

492 

TO 

656 

9 

1.2 

TO 

1.5 

19 

200 

TO 

250 

654 

TO 

820 

4 

1.3 

TO 

1.4 

51 

250 

TO 

300 

820 

TO 

984 

4 

1.4 

TO 

l.s 

14 

300 

TO 

350 

98« 

TO 

1146 

2 

1.5 

TO 

1.6 

20 

350 

TO 

400 

1  144 

Tn 

1312 

l 

1.6 

TO 

1.7 

29 

400 

TO 

SOO 

1312 

TO 

1640 

1 

1.7 

TO 

1.8 

18 

500 

TO 

600 

J  640 

TO 

1908 

2 

1.8 

TO 

1.9 

12 

600 

TO 

700 

1964 

TO 

2296 

2 

1  .9 

TO 

2.0 

8 

700 

TO 

800 

2296 

TO 

2624 

0 

2.0 

TO 

2.2 

12 

800 

TO 

900 

2624 

TO 

2952 

2 

2.2 

TO 

2.“ 

10 

900 

TO 

1000 

295? 

TO 

3280 

1 

2.4 

TO 

2.6 

a 

1000 

TO 

1240 

324  0 

TO 

3937 

0 

2.6 

TO 

2.8 

3 

1200 

TO 

t«no 

3937 

TO 

4593 

0 

2.8 

TO 

3.0 

5 

1400 

TO 

1600 

459  3 

TO 

5249 

0 

3.0 

TO 

3.5 

2 

1600 

TO 

2000 

5244 

TO 

65M 

0 

3.5 

TO 

4.0 

2 

OVER 

2000 

OVFR 

65b  l 

1 

OVCR 

4.0 

3 

o 


11-223 


646 


MILLCREEK 

TEMPERATURE  THRESHOLD  -  2.50CT 
LflMBOfl-  9.3  TO  11.7  wM 


11-224 


MILL  CREEK 


DISTRIRUTTOm  of  RFCOfi*tZ£U  TC“R.  THRESHOLDS  FOR  O  ■  2.5 


5Y  area 


SQUARE 

hfters 

frequency 

0.0 

TO 

100.0 

229 

lon.n 

TO 

20  0.0 

16 

200.0 

TO 

S'*0.0 

11 

500.0 

TU 

1000.0 

2 

lOOQ.O 

TO 

15»0.0 

3 

1500.0 

TU 

2»H'0.0 

2 

2000.0 

TO 

2500.0 

1 

2500.0 

TO 

3oon.O 

1  ‘ 

3000.9 

TO 

one o.O 

0 

oooo.n 

TO 

5000. 0 

0 

5000.0 

TO 

6000. 0 

0 

6000.0 

TU 

8000,0 

0 

8000.0 

TO 

1 0000.0 

0 

10000.0 

TO 

1 5000.0 

1 

15000.0 

TO 

20000.0 

0 

20uoo,o 

TO 

00000. 0 

0 

00000.0 

TO 

8  0  0  0  0 . 0 

0 

80000.  o 

TO 

160000. 0 

0 

OVER 

16000 0.0 

0 

TOTAL  MUMflEU 

UF  TFmP, 

T*Rs  266 

flY  PFRIMtTER  BY  SUAPF 


HCTFRS 

F 

EFT 

FREQUENCY 

SHAPE 

:  FACTOR 

FREQUENCY 

0 

TU 

50 

0 

TO 

16  0 

223 

0.0 

TU 

1.0 

0 

50 

TU 

100 

160 

TO 

328 

22 

1.0 

TO 

1.1 

0 

100 

TO 

150 

328 

TD 

092 

8 

1.1 

Tu 

1.2 

161 

150 

TO 

200 

092 

TU 

656 

3 

1.? 

m 

1.3 

9 

2i>0 

Tn 

250 

656 

TO 
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MILL  CREEK 


#  i 

Va/ 

DISTRIBUTION  OF  Recur,.'!)  ZEO  T£mP.  THRESHOLDS  FOR  O  -  3.0 

BY  AREA 

SQUARE  METFRS  FREQUENCY 
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ARTIFICIAL  COLOR  COMPOSITE  IMAGE  FOR  FLINT-1  AND  THE 
PSEUDO-IMAGE  GENERATED  USING  AREA/ INTENSITY  STATISTICS 


No.  of  Ellipses 
In  Pseudo-Image 


Color 


Intensity  Range 
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Yellow 


Blue 


Below  la 


Brown 


Mean 


